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Abstract 


The assembly of a portable, low-cost, microcomputer- 
based data acquisition system intended for the study of 
atmospheric turbulence is described. The system measures 
fluctuating temperature and wind with a platinum-wire 
thermometer, an inexpensive sonic anemometer and a pair of 
propeller anemometers. The performance of these sensors, 
with particular emphasis on frequency response, is examined. 
Analog conditioning, multiplexing, and 8-bit analog-to- 
digital circuitry, developed to interface the sensors to a 
low-cost RCA VIP-711 training microcomputer, are decribed. 
The microcomputer controls data collection at from 5 to 20 
Hz, buffering up to 3328 bytes. Data collection is stopped 
while data blocks are transfered to standard audio cassette 
tape at 100 bytes per second. The microcomputer also 
retrieves the data from tape and transmits it to a host 
computer at 600 baud via an RS-232C interface. Documented 
assembly and machine language listings of all software are 
included. Suggestions are made for modifying the software to 
allow continuous data collection, variable sampling 
frequencies, faster data transfers to the host computer, and 
real time output. 

The potential of the prototype system is demonstrated 
with examples of atmospheric spectra computed from data 
collected at a moderately uniform and horizontal site. 


Resulting power spectra and cospectra compare well with 
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results from others as well as with theoretical expectations 
in the inertial subrange. The effects of analog to digital 
resolution on the data are also examined. Eight bits are 
sufficient provided analog amplification and offset settings 


are properly used. 
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Chapter 1 
INTRODUCTION 
Our understanding of the atmospheric surface layer has 
been greatly enlarged by the study of careful measurements 
of wind and temperature fluctuations over flat, uniform 
terrain. Parameters of particular interest include the 


vertical flux of heat, 
H= pCpw'd', Ciel) 
and the downwards momentum flux or surface shear stress, 
ace SPEC na Gi) 


where the overbars denote time averages and primes indicate 
fluctuations about a zero mean. In addition to time domain 
Statistics such as these, much emphasis is now placed on the 
Spectra of atmospheric turbulence. There is substantial 
experimental confirmation of the predicted inertial subrange 
behaviour in one-dimensional Eulerian spectra at frequencies 


above n= U/z: 


Shin) =) Oe a eane ea Ty ae, ias:) 
Syini= Swini= (4739Su (0), Gilera’) 
Sah = 0B CU/2ieege We eed eT tras Bree 


The Kolmogoroff constants, @ and £, are relatively well 
established, thus dissipation rates of energy, e, and 
temperature variance, N, may be evaluated from measurements 


made at frequencies far below the actual dissipation range. 
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It is only in recent years that instruments capable of 
making reliable measurements of atmospheric turbulence at 
frequencies in the inertial subrange have become generally 
affordable. For the the all-important vertical wind 
measurements, the sonic anemometer is currently the 
instrument of choice. The basic principle of the device is 
Simple; the difference in times for sound waves to travel in 
opposite directions along a fixed path is proportional to 
the wind speed along the path. However, the first commercial 
sonic anemometers, introduced in the late sixties, were 
complex and costly and therefore not widely used. Recently, 
following the work of Campbell and Unsworth (1979) and 
Larsen at al. (1979), single-component sonic anemometers for 
vertical wind measurements have become available for around 
$1000.' This is only a fraction of the cost of more 
established units and work by Shuttleworth et al. (1982) 
indicates further price reductions and performance 
improvements are still to come. 

Flow ee caused by the acoustic transducers of 
a sonic anemometer become a severe problem when horizontal 
wind measurements are attempted. Consequently, horizontal 
sonics seem likely to remain as very specialized (and 
expensive) research tools. The ionic anemometer, recently 
developed by Barat et al. (1982), may prove to be the 
low-cost answer to measuring horizontal wind fluctuations. 


This simple anemometer, based on the drift of ions created 


‘All prices are Canadian (1982) 
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during a corona discharge, is well suited to both vertical 
and horizontal turbulent wind measurements. At present, 
however, relatively low-cost, commercially-available 
instruments used to obtain horizontal wind variations 
include propeller anemometers (Hicks, 1972) and drag 
anemometers (Smith, 1980). 

Much of the research aimed at developing low-cost 
turbulence instruments centers around an interest in 
obtaining reliable, routine (possibly even remote) flux 
measurements by use of (1.1) and (1.2). This so-called eddy 
correlation method has attracted much interest from workers 
in agriculture and forestry, as well as meteorology, who 
have traditionally used highly empirical, aerodynamic 
methods to estimate fluxes from profile measurements. This 
emphasis on flux meaSurements is reflected in the data 
acquisition products currently on the market. Still 
following an approach first used by Dyer et al. (1967), 
commercial units tend to be dedicated to computing and 
displaying real-time covariances. While analog circuitry was 
originally used for this purpose, modern systems are 
microprocessor-based and perform digital calculations. 
Either way, the drawback with such an approach is that it 
precludessany further analysis ofethe data. As a result; the 
study of atmospheric spectra continues to be done with 
on-site minicomputers (e.g. Redford et al., 1981) or 
multichannel FM analog recorders which are later digitized 


with a minicomputer-based facility. Both approaches are 
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expensive, cumbersome, and require a high degree of 
technical expertise. 

For conventional meteorological applications, micro- 
processor-based data logging is becoming commonplace. The 
least expensive systems, which store data on audio cassette 
tape, start at about $2000. Retrieving and transmitting the 
data to a host computer for processing requires another 
device of comparable cost. Weihofen and Woehl (1981) have 
shown, however, that a commercial training microcomputer, 
costing less than $200, could be used in place of such a 
system. Their device saves data on audio cassette tape using 
the 300 baud (30 bytes per second) Kansas City Standard 
format. By using this popular audio encoding scheme, the 
tapes of Weihofen and Woehl could be read directly by the 
host computer system available for their data processing. 
Unfortunately, use of the Kansas City Standard results in a 
system that is too slow for collecting multiple channels of 
turbulence data. 

The slow speed of the Kansas City format has led many 
microcomputer manufacturers to abandon it in favour of their 
own faster versions. Several training computers, in the 
class used by Weihofen and Woehl, are currently available 
with cassette interfaces operating at around 100 bytes per 
second which makes them suitable for collecting turbulence 
data. Such non-standard formats, however, are unlikely be 
supported by the host computer to be used for data 


processing. This means the microcomputer must be used to 
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control data retrieval as well as data collection. 

The objective of this project was to assemble a 
microcomputer-based data acquisition system capable of 
collecting several hours of turbulence data and transferring 
these data to a host computer for processing. In addition to 
being inexpensive, it was desired that the system be easy to 
construct without the need for specialized development 
equipment. In order to accomplish this the system has been 
kept as simple as possible. The possibility of combining 
real-time calculations with the primary task of data 


collection, although desirable, has not been explored here. 
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Chapter 2 


A LOW-COST INSTRUMENT ARRAY 


2.1 Instrument Requirements 

SENSOrsS are the Starting point of any data acquisition: 
system. Those used to measure atmospheric turbulence require 
much better resolution and frequency response than are 
needed in conventional meteorological applications. As well, 
turbulence sensors must be carefully designed so that they 
ao not Significantly modify the characteristics of the flow 
which they are attempting to measure. 

The response of an instrument to a Sinusoidal input may 
be characterized by an amplitude function, H(n), and a phase 
faq, ¢°n),ssuch that sa sinusoidal input’ signal, sin(27nt), 
is transformed to H(n)sin{2mnt-¢(n)}. Here, it is more 


convenient to use the transfer function 
Tin)= He Gnp= S' (n)7S(n), (2-34) 


where S'(n) is the computed estimate of the actual power 
spectrum, S(n). Knowledge of T(n) reveals the limitations of 
an instrument and, up to a certain extent, enables them to 
be corrected. Typically, the half-power point, where 
T(n)=1/2, is taken as the representative cut-off frequency 
of a sensor. Figure 2.1, from McBean (1972), permits one to 
estimate the loss of vertical flux due to a high frequency 
cut-off as a function of height, wind speed, and stability. 
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Figures2ak..sApproximate loss of vertical flux due toa 
sharp, high frequency cut-off (after McBean (1972)). 
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Surface require cut-offs ranging from 1 to 5 Hz in order to 
keep the loss of flux under 10%. 

Instrumental phase response is often ignored because 
power spectra are independent of it. Differences in ¢(n) 
between two instruments will, however, effect covariances 
and cospectra. A complex cross-spectrum, S;2, may be 
corrected for instrumental response by subtracting Ag=¢2-¢, 
from the phase and dividing the magnitude by by H;,H2. The 
cospectrum may then be obtained from Re{S,2} and numerically 
integrated to give an improved estimate of the covariance. 

In the present project two propeller anemometers, a 
sonic anemometer, and a platinum wire thermometer are used 
to measure atmospheric turbulence. Each instrument 
represents a different trade-off between cost and 


performance which must be individually examined. 


2.2 The Propeller Anemometer 

The propeller anemometer, developed in the early 
sixties, is a durable, low-cost, fairly sensitive wind 
sensor. The widely-used Gill UVW array (consisting of three 
propeller anemometers mounted at right angles to each other) 
presently sells for around $1500.? Each anemometer consists 
6flandightweight, *four=blade;,*helicoid propeller? (23 cm 
diameter) directly driving a d.c. generator which produces 
approximately 50 mV/(m/s). Ideally, a fixed propeller 


anemometer would respond only to the component of the wind 


2Made by R. M. Young Co., Traverse City, Michigan 
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vector parallel to its axis, thus allowing a fixed array of 
three orthogonal propellers to completely determine the wind 
field. In practice, Gill propeller anemometers deviate from 
this so-called cosine response as the angle of attack (angle 
between the rotation axis and the wind vector) increases. 
Drinkrow (1972) and Horst (1973a) have investigated this 
problem and devised iterative schemes for recovering the 
true wind components from those meaSured by an orthogonal 
array of propellers. 

The frequency response of the Gill anemometer is of 
more serious concern. Typically, it is modelled with a 


transfer function of the form 


To(ni=l7 tte Ceanb/ Us) (2,2) 
and a phase lag 


?o(n)= tan” '(2mnL/U), b233:) 


where L is called the length constant. By examining the 
response of propellers released from rest in steady, wind 


tunnel flows, Hicks (1972) suggested that 
L= Becosi. +o, C245) 


where Lo is the length constant measured with the mean wind 
aligned with the anemometer (about 1.0 m for a 23 cm 
propeller) and © is the angle of attack. A half-power point 
at n= U/2mL means that this sensor will frequently lose 20% 
or more of the vertical flux and rarely allow dissipation 


estimates from the inertial subrange without large 
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corrections to computed spectra. Figure 2.2 shows the large 
variations in T (n) due to changes in axial wind speed. In 
addition, individual variations due to bearing wear and dirt 
can also substantially degrade the frequency response. This 
makes applying meaningful corrections to propeller spectra 
difficult, since response lengths based on individual 
calibrations should be used, rather than those supplied by 
the manufacturer. 

If one propeller is aligned directly with the mean wind 
(O©=0°) for optimum performance, (2.3) implies that the 
frequency response of an orthogonal propeller (lateral or 
vertical) will be extremely poor. Fortunately, Fichtl and 
Kumar (1974) have shown that the situation is not quite so 
bad in a turbulent atmosphere. For angles of attack near 90° 
they found that the response is primarily a function of the 
intensity of the turbulence along the propeller axis, 
resulting in length constants typically 4 or 5 times less 
than those given by (2.4). Even so, the frequency response 
is continuously varying and the transverse propeller is 
almost always operating in its non-linear regime below 
1 m/s. Stalling occurs frequently as the axial wind speed 
falls below 20 cm/s. 

For horizontal wind measurements, mounting a perpen- 
dicular parr of propellers such that theymean horizontal 
wind makes an angle of 45° with each, seems to be the best 
approach: In thiswspecialicasey bothepropeilersahave: the 


same response (L= 1.2°m) and a single transfer function can 
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Figure 2.2. Transfer function (a) and phase lag (b) of a 
Gill propeller anemometer with L=1 m. The abscissa may also 
be read as reduced wavenumber, f, in which case the wind 
speed labels may be read as heights (i.e. 1 and 5m). 
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be used to modify the Spectra of the resolved longitudinal 
and Deeeeal components. All other cases will involve 
combining data from sensors with different responses. Brooks 
(1977) derived an effective time constant, using (2.4), for 
a wind component resolved from a pair of perpendicular 


propeller anemometers, given by 
L'= L{1.093+.093sin(40-7/2)}. C226) 


Despite the limitations of the propeller anemometer, 
low-cost and durability have resulted in its widespread use. 
The lack of an inexpensive alternative has resulted in the 


use of a pair of horizontal Gill anemometers here. 


2.3 A Low-Cost Sonic Anemometer 

One-dimensional, sonic anemometers are well suited to 
vertical wind measurements because of their relatively good 
frequency response, linearity, and absence of a threshold. 
Flow distortions in the wakes of transducers, which greatly 
restrict the allowable angle of attack of horizontal sonic 
anemometers, are not a problem for vertically-mounted units 
Since the mean wind is almost always transverse to the sound 
paths The CA27)> a commercial Wwersiom jot Aphe continucus-wave 
sonic anemometer designed by Campbell and Unsworth (1979), 
was purchased for this project, at a cost of about $1200. 
The instrument is shown in Figure 2.3. Most previous sonics 


have used separate transmitters and receivers at each end of 


>Made by Campbell Scientific, Logan, Utah 
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Figure 2.3. The CA27 Sonic Anemometer. Only a single ultra- 
sonic transducer 1s used at each end of the 10 cm sound 
Patni. 
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the sound path. In the CA27, piezoelectric transducers 
interchange the tasks of transmitting and receiving ultra- 
sound a hundred times per second. The unit operates from a 
12 V battery, has an output of 1 V/(m/s) and a range of 
+3 m/s. The effect of temperature variations on both the 
speed of sound and the characteristics of the transducers is 
Mot caxen into account .by the” circuitry of this sonic 
anemometer, thus limiting accuracy to about +1 cm/s. 
Although susceptible to rain and dew, the small size of 
the transducers employed by the CA27 allow the use of a 
sound path only 10 cm long, half the distance used by most 
sonics. This is significant, as it is the averaging of the 
wind over this path which ultimately limits the high 
frequency response of the instrument. For the case of mean 
wind along a sound path of length d, Mitsuta (1961) derived 


the transfer function 
T.,(nJ= ((U/and)sintmnda/u) } 7 (255) 


by applying Taylor's frozen turbulence hypothesis, x=Ut, and 
integrating a sinusoidal wave over the time interval 
appropriate to the sound path. This approach also indicates 
that there is no phase shift associated with the horizontal 
line averaging process. A comparable analytic expression for 
a vertically mounted sonic anemometer has not been found, 
although Kaimal et al. (1968) and Horst (1973b) have 
produced numerical approximations based on the assumption of 


an inertial subrange. The following fifth-order polynomial 
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closely fits their tabulated data for a vertical sonic 


anemometer with a single, 10 cm sound path: ‘ 


Ts3(n)= exp{-0.0268-0.0664x-0.0662x?-0.0201x3 
=0.000269x°+0 .0000383x°%, C276) 


where x=ln(n/U). Plots of Ts; are shown in Figure 2.4. 
Although much better than that of a Gill propeller, the 
frequency response of the sonic anemometer still has a 
strong dependence on the mean wind. However, this correction 
is relatively small in the range of interest here and, since 
U is the only important meteorological variable involved, 
Ts3; can be applied to spectral data with a high degree of 


confidence. 


2.4 A Platinum Wire Thermometer 

Fluctuations of temperature are much easier to measure 
than those of wind. Generally, thermocouples or resistance 
wires with small size and mass are used. A fine-wire 
thermocouple (Campbell Scientific 127), intended for use 
with the CA27 sonic anemometer, sells for over $300. Here, 
following Haugen et al. (1967), a single strand of 99.9999% 
pure platinum wire, 25 um in diameter, was loosely strung 
parallel to the sonics path. This sensor, and accompanying 
GUreusrenypecost about $10.5 


‘TheNes values, (labeledwts)ipin, Table 3 of Kaimalyet al. 
(1968) should each be displaced to the next-lowest tabulated 
freotency., Thergrape Of Topigerigure J Ol that paper. iS, 
however, correct. 

5Platinum wire was obtained from Alpha Co., New York. Cost 
(imeearly J9G62),.Was $3/ perlumerer, 
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SONIC TRANSFER FUNCTION, Ts 
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Fioure 2.4. Transfer Ttunct yonvot a .verticalysonic ranemo— 
meter with a 10 cm sound path. Also shown (dashed line) is 
theltranster suinction Lomrasnonmizontal sonic with a 10.cm 
Sound patina below wo whz the two-eurvesgqenerally -agree*quite 


well. 
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The advantages of platinum are well known and include 
good stability and linearity over temperatures of 
meteorological interest. The change in resistance with 
temperature, 0.00395 R(0°C) per °C, is large compared that 
of most metals. The resistivity of pure platinum is 

10.6 wohm-cm which, for the chosen diameter of 25 um, gives 
a resistance of 2.16 ohm/cm. Kaimal (1968) has shown that 
the frequency response of such thin wire is predominantely 
determined by the line averaging process; therefore, the 
transfer function presented for the vertical sonic 
anemometer is directly applicable. 

A 10 cm strand of the platinum wire (21.6 ohms) was 
placed in the Wheatstone bridge shown in Figure 2.5. A ten- 
turn, 10 ohm, variable resistor, Ra, iS adjusted to bring 
the resistance of the platinum wire plus leads up to 28.0 
ohms. This is easily done by nulling the bridge output, eo. 
The fixed resistors (metal-film type with low temperature 
coefficient) were matched to better than 0.5%. Analysis of 
the’ circuit in Figure 245, under the assumption that AR, the 
variation of the platinum wire resistance, is much less than 


R, gives 
@o= E,AR/4R C287) 


In order to avoid self-heating and wind speed dependence, 
the current through the platinum wire must be kept low. 
Using B,=5 V and R,=1 kilohm gives e,=135 mV and a current 


of. 2.4.mA through the.wire. Power dissipated by the wire is 
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Figure 2.5. The platinum wire thermometer. The resistance 
Of a 10 em strand of 25 um platinum wire,2s :about? 22)-ohms. 
The temperature bridge adjust, Ra, 1S used to bring the 
active arm up to 28 ohms in order to balance the bridge. 
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approximately 150 uwW. E, is provided by a 7805 voltage 
Feguiator' which varies by less’ than 1% from ‘0°C to 30°C. 
The bridge output, eo, is very small, and must be 
applied to a high-gain, low-noise, differential amplifier. 
Using a gain of 5000, a 22 ohm strand of platinum wire was 
calibrated in a water bath which gives more stable and 
easily controlled temperatures than air. Using the range 20 
£07235; cF) 6 icalibration coefficient fof] 0.107+0+004 mVZ°C was 
obtained. The theoretical value for a 22 ohm wire is 
0.105 mv/°C. From this good agreement it appears that a 
Simple calculation based on the resistance of the platinum 
wire is sufficient for calibration. The calibration 
coefficient, in mV/°C, is obtained by multiplying the 
resistance of the platinum wire at room temperature (in 
ohms) by 0.00476. Figure 2.6, a comparison between the 
platinum wire thermometer and a commercial thermocouple 
(Campbell Scientific 127), clearly shows a high visual 
correlation. The effect of line averaging can be seen in the 


Slightly smoother platinum wire trace. 


2.5 Instrument Mounting 

Mounting the selected instruments in the field is not a 
trivial task. The continuing controversy over the results of 
the Kansas experiment (e.g. Wyngaard et al., 1982) 
underscores the importance of avoiding obstructions which 
may distort the flow being measured. Proper alignment of the 


wind sensors is also vitally important. Wieringa (1972) 


<! 


= 


apaatov. 20Bt 6 ae sts sihcaed a jee 08% 
MO mals, qr ede awa yd-aeirey’ dobde'® 


OE oF oh 
ad. 2G 


his ! Tfens viet ei ae wetuc spbiad ~~ 


oki 


ve 


eA atte 


_eslttlame leftraasit is. (se fon=wal pap dpidos oF 4 
25W 4 mintialqa to hfetde ade. Gee OO0e: Lo riled <a 
:[cbt% atom aevig bide deed >egew a of Bese 
apres si Ake 2a leet sag aS iuisis cus? bs epiaitiati a) 
@ Vin. £00,028 0 : HS Lollies. sid gee Ewe & i 
axtv to $8 5) 791 au tee lee é7eaowds i 
E 2602 Ayseqqe 22 itesertee BoeP pity ore. 
wwiTelq sds Io-saneseies7 ay? We beac a: teluoteae 
ae 
. gintias sit) , re iserdel ss Tee siete et 3ue ate = 
sis oniviqrsaon va Ger erat Se 2° we at 00} oH a 
Siu jis Tegnss fice: ts) sateen iseag ete Te sonet2i aii 
F Vi 
3. dueweed iNeet y6cme> &.)a.5 enyor’ Be ¢ \e tm nike 
2 we DIM Srl tee oTenhee! s hive Jessmotadits' sities mun 
toei+) eid 6 3wone Vliyaeio WOES sikiSasind “Ie ams); 
en4 A 4 80 eo Bil peteve Sarr ‘tg time edt colts teres 
2° SBShee fi0W nuintte der staan ulaigily 
; : be 
Shad hice Pertatt ri ae 
’ 
en: et, mig? gaz ei. saris 1724t Gerssal<a ‘en: onbtagom : 
baie 


ESiuee, sf? evo Yetars 23G6 


{R007 |, Lage Siseon yy. 2p.) ‘deqthegns conn 


ait ‘peottowagede’ pribievs 


dt to: inatings is" wagons aBetianen orhee, wert i ste aeee 
eheete ene. basi oats at ax ude 


| at 


ae 
te S 
— "~~ 
aa oP 2 . art o 7 
_ a le 7 a 
a ~. i 


I-EMUNISioS. sit .deee ising 


io Sugeno ets ai} cesesetetma ' 
7 


] Py f 


a 


a” 


20 


Thermocouple 


HT I/| | OSC 


| 10 seconds | 


Platinum Wire i 


0.552: 


Figure 2.6. Comparison between platinum wire thermometer 
and a commercial, fine-wire thermocouple. The traces were 
obtained in the laboratory by gently blowing across the 
probes. The thermocouple was located about 2 mm from the 
center of the resistance wire. A high viSual correlation is 
evident although the thermocouple seems to show more fine 
detail. 
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found that one degree of sensor tilt could result in errors 
of up to 10% in momentum flux and 5% in heat flux. 

The CA27 sonic anemometer seems to be poorly designed 
in terms of minimizing obstruction to the flow. Hogstrom 
(1982) has recently found errors of 10 to 20% in variances 
and covariances caused by flow distortion around a vane- 
mounted hot-wire turbulence sensor which is much more 
aerodynamically shaped that the sonic anemometer used here. 
It therefore seems likely that distortion of the flow due to 
the cylindrical base of the CA27 could significantly effect 
the vertical wind measurement. To reduce this problem, the 
sonic probe was remounted on a seven pin, MS-type, Amphenol 
connector compatible with the standard Gill UVW array. The 
filter and blower housed in the base of the Gill array (used 
only for continuous operation applications with AC power) 
were replaced with the electronics of the sonic anemometer. 
The platinum wire probe was mounted on the sonic probe 
parallel to the sound path. The resulting wind-temperature 
array, shown in Figure 2.7, reduces obstruction problems 
while providing a single unit which is fairly easy to handle 
and level. Leveling by spirit levels and visual inspection 
Canetypically be done tojwithin ]ti)-. When sonly vaysingle 
level of measurement is attempted, the array should be 


mounted as high as possible in order to optimize sensor 


response. 
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Figure 2.7. The wind-temperature array. The sonic 
anemometer and platinum wire thermometer mount into the 
socket ‘on the Gill UVW array originally intended for a 
vertical propeller. 
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Chapter 3 


ANALOG SIGNAL CONDITIONING 


3.1 The role of Signal Conditioning 

The outputs from the instruments, eventually to be 
digitized, must first undergo considerable analog 
processing. Typically, the small fluctuations of interest 
will be superimposed on a large d.c. component representing 
the mean value of the signal. This must be removed and the 
fluctuations must then be amplified so that they utilize the 
available digitizing range (nominally 0 to 5 V for all 
channels) as fully as possible. High-pass filters can be 
used to remove d.c. and low-frequency components of the 
Signals, although, a simple offset, if possible, is 
preferred as it allows easy reconstruction .of the absolute 
Signal level. Finally, it 1s important that the signals 
undergo high-pass filtering in order to reduce errors due to 


aliasing. 


3.2 Aliasing 

Aliasing errors are introduced by the digitizing 
process. Only frequencies less than half the sampling 
frequency, the so called Nyquist frequency, no= 1/2At, can 
be resolved. Energy from the higher frequencies is folded 


back around the Nyquist frequency according to 


Sa(n)= Se nlee S(2mno-n)+S(2mno+n), (aa) 
m=1 
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where Sq the aliased estimate of the true power spectrum, S. 
Applying a rectangular filter prior to digitizing, which 
removes all frequencies above no, would eliminate aliasing 
completely, but such ideal filters are hard to come by. A 
practical, two-pole, low-pass, Butterworth filter with 


half-power point at no has a transfer function given by 


Beene =1/1 +n ne (3.2) 
and phase 


o;(n)= tan \inen/(ne4an*)}. (3.3) 


These are shown in Figure 3.1. The phase lag of the filter 
can be ignored only if all of the channels experience the 
same filtering. 

Given the natural n°’? roll-off in the inertial 
subrange of most atmospheric spectra, a two-pole filter, 
although unable to completely eliminate aliasing, will 
confine the contamination to a single fold (negligible 


Céntuibution from above 2n,). 
Soin] 2 (nS) teen )/S tana). (3.4) 


As can be seen in Figure 3.2, S;(n) is a good approximation 
to the true spectra due to the close cancellation of 
abvasing by the two-pole filter .®or further improve the 
spectral estimate it is necessary to provide an explicit 
expression for S$(2n.-n). ‘This can be done by extrapolating 
the measured spectra (not including the aliased end) or, 


followinguhkaimal et.al. (i772) ,sinertial” subrange behaviour 
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LOW-PASS TRANSFER FUNCTION 
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POWER SPECTRA (Arbitrary Units) 
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Figure 3.2.” Btfect of a two-pole, low-pass filter’ on 
aliasing in the inertial subrange. Sa represents the aliased 
estimate of S in the absence of any filtering. S; Shows the 
result of low-pass filtering with half-power point at the 
Nyquist frequency (5 Hz). Finally, Sc shows the result of 
SCOrlechiicGm ine: filtered Soecirogsio . OD UCne ser rec. s ot 
low-pass filtering. While Sc is clearly an improvement over 
Sa, the best overall estimate appears to be Sy. 
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could be assumed. In most cases, however, this extra work 


will not appreciably improve the spectra. 


270 Practical Circuitry 

To provide the required amplification, offset, and 
filtering,the output from each instrument was passed through 
Chem precondition circuit shownvin Figurev3é3. It consists of 
a differential preamplifier followed by a high-pass filter, 
an inverting amplifier with adjustable gain and offset, and 
a low-pass filter. 

The major function of the differential amplifier is to 
reference all the incoming signals to a common ground. By 
reversing the role of the input terminals, this stage can 
also be used to set the final polarity. For most sensors 
R,= R2= R3= Ry= 1 megohm, thus giving unit gain and input 
impedance of about 500 kilohm (all resistors were matched to 
better than 0.5%). For the temperature signal, R, and R; 
were replaced with 100 kilohm resistors in order to give 
this stage a gain of 100. The resulting input impedance of 
only about 50 kilohms is not a problem, due tothe low 
output impedance of the temperature bridge (less than 30 
ohms). 

The output from the differential amplifier then goes 
into a two-pole, Butterworth, migh-pass filter with a halt- 
power point at 0.005 Hz. The ‘output from this filtes 
connects into the inverting amplifier stage via the toggle 
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Differential 
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amplifier R13 R14 
Offset adjust MS es pee 
R9 Ply 10k R10 
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° low-pass Low-pass filter 
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Option 5 Hz 
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200k 2k 2k 
10k 10k 10k 


10k 10k 10k 

20k 5k 5k 

10k 10k 10k Spare 
10k 8.5k 8.5k j|channels 


Figure 3.3. The prototype analog precondition circuit. 
Component values which vary from channel to channel are 
shown in the table. All resistors are 1% metal-film types. 
The operational amplifiers are each 1/4 of a TLO84CN. 
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of the circuit if desired, in which case the output from the 
differential amplifier is sent directly to the inverting 
amplifier. The switch does not close off the input to the 
high-pass filter as this would require the filter to 
resettle (approximately 200 seconds) whenever S, is thrown. 

The gain of the inverting amplifier stage is given by 
Re/R7 where Re, constructed from a 10-position rotary switch 
and 100-kilohm resistors (1%), is variable from 100 kilohm 
to 1 megohm in 100-kilohm steps. The value of R; depends on 
the sensor. For the sonic anemometer R,= 200 kilohm has been 
Pouna-to- work well;Yqiving’ gains’ from: 075”to"5; “For 
propeller anemometers R7= 20 kilohm results in acceptable 
gains from 5 to 50. The platinum wire thermometer amplifier 
uses R,= 2 kilohm, which, in conjuction with the gain of 100 
in the differential amplifier, gives gains from 5000 to 
50,000. The ranges were selected, on the basis of field 
tests, so that the typical range of turbulence levels can be 
accommodated with reasonable sensitivity. 

The adjustable potentiometer arrangement attached to 
the noninverting terminal of the inverting amplifier stage 
allows a d.c. offset to be added to the signal, in order to 
approximately center it in the 0 to 5 V digitizing range. A 
ten-turn potentiometer was used for P, to allow good control 
over a wide range. 

Finally, the signal is subjected to low-pass filtering. 
A two-pole, Butterworth filter with half-power point at 5 Hz 


is appropriate for the typical sampling frequency of 10 Hz. 
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However, to accommodate different sampling rates, the switch 
S2 provides a path to an optional, interchangeable filter. 
In the prototype system, the circuit shown in 
Figure 3.3, excluding the optional filter, was implemented 
On Standard, one-sided, 22-pin, printed-circuit cards, each 
using a single, inexpensive, monolithic quad operational 
amplifier (TLO84CN). Temperature dependent offset drifts are 
not critical in this application as the analysis of data 
invariably requires the removal of mean values and linear 
trends over a few minutes. The cost of parts, per channel, 
was about $15 although this could vary substantially 
depending on the quality of the hardware selected. 
Commercial quality integrated circuits, which have a 
recommended operating range of 0°C to 70°C, were generally 
used here. The prototype system includes two spare channels 
for a total of six. Optional filters for all six channels 
were implemented on a Single printed-circuit card. The 
operational amplifiers in the system are powered by +8 V 
obtained from a pair of 12 V lantern batteries via common 


monolithic voltage regulators (LM317,LM337). 


3.4 Field Operation 

The six channel prototype Analog Precondition Unit 
(APU) is shown in Figure 3.5. Two 12 V lantern batteries, 
which supply the power to the unit and the sonic anemometer, 
fit inside the unit. The basic wind-temperature array plugs 


into the rear panel of the APU with a single connector. 
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Separate input connectors are provided for the two spare 
channels, X and Y. The channel switch allows any of the 
conditioned signals to be displayed on the analog meter. The 
offset switch allows reading the applied d.c. offsets. 

In the field, it is necessary to adjust the gain and 
offset settings of each channel to obtain signals with means 
of about 2.5 V (25 wA on the analog meter) and fluctuations 
between 0 and 5 V (0 and 50 wA). The analog meter can be 
used for this purpose, however, a sStrip-chart recorder works 
pivchebetter, Thesrecorder,»should be.set»on a.0.to.5.V range 
and connected to the meter output at the rear of the unit. 
If additional strip-chart recorders are available, 
individual channels can be monitored via conditioned signal 
outputs provided at the rear of the APU. 

During field testing, it was found that high-pass 
filtering was never required for vertical wind measurements. 
For temperature, a fixed offset could generally be used for 
intervals of only about 15 minutes. The best gain settings 
at the start of an experiment were such that the maximum 
fluctuations covered from 50 to 75% of the available range. 
A slightly larger gain could be used if high-pass filtering 
was employed; however, filtering was avoided whenever 
possible due to the limitations it imposes on the low 
frequency contributions to fluxes and spectra. In addition, 
high-pass filter corrections are somewhat uncertain due to 


the difficulty of trimming and testing such low frequency 
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Chapter 4 


DIGITAL DATA COLLECTION 


4.1 A Low-Cost Microcomputer 

Following the most inexpensive approach possible, the 
conditioned analog signals were channeled through a multi- 
plexer (MUX) into a single analog-to-digital converter 
(ADC). In order to operate correctly, the MUX and the ADC 
have to be provided with various control signals. The 
digital data produced by the ADC must then be collected, 
encoded as audio signals, and recorded on a cassette tape. 
These tasks may be performed easily and economically with 
the aid of a commercial microcomputer. The approach used in 
the present system is to collect the data temporarily ina 
reserved memory area in the computer. When this buffer is 
filled the collection process is stopped and the data are 
transferred to a cassette recorder. Data collection resumes 
when the tape transfer 1s complete and the cycle continues 
for as long as required. This is clearly not the best 
approach; continuous data collection would be preferable. 
However, the approach adopted is the simplest one possible 
and, therefore, was considered an appropriate starting 
point. 

Numerous low-cost, single-board microcomputers, 
intended primarily for educational and recreational use, 
were available in early 1982 for under $300. The RCA COSMAC 


VIP™CDPI8S711, Shown in Figure 4.1, was chosen for this 
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project although, unfortunately, it has since been 
discontinued by RCA. It is based on the 8-bit, COSMAC 1802 
microprocessor running at 1.7609 MHz. The CMOS (Complemen- 
tary Metal Oxide Semiconductor) fabrication of the 1802 and 
other components in the VIP-711 provide low power 
requirements and high noise immunity. On-board memory 
consists of two kilobytes of static RAM (Random Access 
Memory), empty sockets for an additional two kilobytes of 
RAM, and 512 bytes of ROM (Read Only Memory) containing an 
operating system. The VIP-711 also includes a cassette 
interface, video interface (an RF modulator is needed to use 
a TV without a video input), parallel 1/0 port, expansion 
connector, hexadecimal keypad, tone circuit, and a power 
Supply. Except for the power supply, the system fits onto a 
Single board measuring about 20 by 30 cm. 

The basic VIP-711 cost $260. Adding a cassette tape 
recorder,* an RF modulator,’ and filling the empty RAM 
sockets brought the price of the useable system to about 
$400. Various accessories which fit into the expansion 
connector were also purchased including a four kilobyte RAM 
expansion board ($150), an EPROM (Eraseable Programable Read 
Only Memory) Stream ($150), and EPROM interface board, 
($60). These relatively expensive options were not needed 
forthescurrent-project. 

The VIP-711's built in cassette interface, implemented 
largely in the operating system software, is vital to the 


Sony TCM=ils ty) OU. 
’"M & R Enterprises Sup "R" Mod-II, $40. 


ny 


need oocke. ead +h a ee : 
Sas: SAMROD .tBfee we: aig: & 2: 7 ) 
-nemsiqnooh Séualcaelt. stake eet 4B, eiomey 
reNag wos ab tet TH eee a igs onan a 
{ronan | Lysod <0 Ni iin BE bSn) dpi bas & ae mM 
2RE2TA OGD AAT MAR Stage Se aedysol ey ‘ows: he al 


Sia S061 $49 fo nofteotadad ere 


ae 


bo aszectol To) ow £1 Siea a Lit SSR Tae tama a 

is pHeh etiod .< Viagem 0 ae beet) UDA f Be ge ye ne. 
t2ie2geo s aebulsrnt obbe Hipbany acer. eure satel 

a2b. af Heese jet. sosslybant ag a) sealaqyed qebiv a4 
tae OAD ‘oliaiey. , Ar teak wabev a rwodstyg 
FovCd) 4 Pas FA o shod seated, Lenk sebexat 7928 “a 
C= SUS 4 ¥eGuie bwag auty 163 eines i 
nial A Te 4: Vode oniieienit a 


eed) Sov bObng 46 BELEBA /.DSES suas. in 
: rer Y BPMLESG | + ZAG ‘hime ied at 


<_< 


ne pataiwes eid aan biame fords hip suai 


MAM. siydel 24) Sy5d te Party bsrettny Calin: 2180 O82 : 
bk on sidemeugert etdebansa FSORTSS es: bere) faced red 


Prbet siathaunt von pre hpete) Sia apotre ct 
keSgan tor arsv anee a s ov ierones eisvitais: est? 


A a 


es 
Seager a7a5 £) 44432 rr py: " 
5 saith Fass ne thas e for~tiy: Pn ents 


oy od 


sift. 32 baghe =! ye vee ein aiid AE a 


si-daacelll 


ha 
a 
he 
ey 
aj 
% 
a 
it 
tl 
‘a 
bee 
~ 
& 
oA 
a 
uv 
Ss 
Co 
eis 
v 
ip 


35 


a 


Video Out 


Power “(5 Volts) 
BhasStic Gover 


b 


vat: EXPANSION CONNECTOR I/O PORT 


0 +VDC 
ae VIDEO 
VIDEO INTEE INTERFACE eer eGo 
Rae I/O INTERFACE aMIC 
VOLT 
OEAR 
ine aos 
eee ileal ie tale Gal 
4K RAM dal eld AND RevEORED We INTERFACE 
beandR[ =e 1802 CPU 6 
Nora, frame STATUS 
© 
beige feed [seme aka HEXADECIMAL 
KEYPAD 
TONE CIRCUIT EXPANSION 
INTERFACE RUN 


SWITCH 


Figure 4.1. The RCA COSMAC VIP-711 microcomputer (a) and 
layout of its major components (b). The VIP-711 comes with a 
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current application. The interface can transfer N pages (256 
bytes per page) of data between RAM and tape in 5 + 2.5 N 
seconds which is almost 100 bytes per second. The first five 
seconds are used to record a constant tone which signals 
Start of a data block during tape playback: This interface, 
while just fast enough for the current application, is over 
three times faster than the Kansas City Standard interface 
used by many single board and home microcomputers. To obtain 
a significantly higher transfer rate than that offered by 
the VIP-711 would require an expensive tape transport or 
digital recorder. 

The VIP-711 must be programmed in COSMAC machine 
language via the hexadecimal keypad. A detailed discussion 
of this language is given in the RCA User Manual for the 
CDP1802 Microprocessor.* The COSMAC instruction set is 
efficient and relatively easy to usSe. There are only 91 
instructions compared, for example, to 153 offered by the 
popular 2-80 microprocessor from Zilog. Although poorly 
suited to mathematical applications, the COSMAC is well 
Suited to the present control task. The operating system and 
video interface, provided with the VIP-711, aid in entering, 
checkingtand = running*programS SinRAM; «unfortunately Prtnere 
are few provisions to aid in program development and 
debugging. Such basic editing features as insertion and 
deletion are absent. Also missing is the option of single 
step program execution offered by some single board systems. 


®Available from.RCA Solid State Division, Somerville, New 
Jersey 
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Consequently, development of software is considerably 
Simplified by access to a COSMAC cross-assembler. This is a 
program, on a host computer, which will accept COSMAC 
assembler language mnemonics’ and translate them into 
machine language hexadecimal code. A limited version of the 
RCA COSMAC cross-assembler is supported by the Computing 
Services at the University of Alberta and was used in this 
project. 

Programming the VIP-711 requires an understanding of 
the register structure of the 1802 which is depicted in 
Figure 4.2. The 8-bit data (D) register, which most 
manufacturers refer to as the accumulator, is used to hold 
the input and output of all arithmetic operations. The 
Single-bit data flag (DF) associated with D indicates 
carries or borrows. The 1802 has sixteen multipurpose, 
16-bit registers not found on most microprocessors. Rather 
than having a special register to act as the program counter 
(which always contains the address of the next instruction 
to be executed), one of the general purpose registers is 
designated program counter by the 4-bit program counter 
pointer, P. Similarily, the data register pointer, X, is 
used to point to a register which will contain the address 
of the next data byte. The 8-bit T register is used to 
temporarily store the current P and X values when an 


interrupt occurs. There are four input or external flags 


"For example, SEQ is used instead of #73 to represent the 
command to set Q on. The # iS uSed to indicate a hexadecimal 


number. 
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Carry Accumulator 


eC) e[L_@ 


Register Matrix 
Sixteen 16-bit registers 


Pointer to 
Program Counter 


Bale (Adeel iat) 


Pointer to 
Data Register 


ta ae 


Temporary storage for X, P 
following an Interrupt 


fips were ee 9 


Interrupt Enable 


aoe 


Upper byte Lower byte 


Internal Flag and 
External Flags 1-bit Output 


EFi[ | EF2( | EF3[ | EF4[ | o{ | 


Figure 4.2. The registers in the COSMAC 1802 microprocessor 
after Peatman (1977). 
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(EF1 to EF4) to provide information from the outside. 
Several branching instructions are associated with each of 
these. There is also a Single output bit, designated by Q, 
which can be turned on or off in software. In the VIP-711, Q 
1s connected to a speaker via a tone circuit. Q is also 
available, in buffered form, on the I/O port. 

The I/O port of the VIP-711 is an important feature not 
offered on all single-board microcomputers. It includes an 
8-bit input port, an 8-bit output port, power connections, 
an input strobe, plus buffered connections to Q, EF3, and 
EF4, The input port allows the #6B instruction to read a 
byte (false= 0 V, true= 5 V) into the data register and the 
memory location designated by X. An external signal applied 
to the input strobe (INST) can automatically cause the byte 
on the input port to be latched. A #68 instruction can then 
read this byte into memory. The #63 instruction places the 
contents of D on the parallel output port, OUTO to OUT7, in 
order to control external equipment. The pin designations of 


the port are given in Table 4.1. 


TABLE 4.1 VIP-711 1/0 PORT CONNECTIONS 


PIN SIGNAL PIN SIGNAL PIN SIGNAL 
A INO (LSB) K INST U OUT6 

B IN1 LE EF4 V OUT7 

C IN2 M OUTO (LSB) Ww Q 

D IN3 N OUT 1 xX EF3 

E IN4 P OUT2 ve +5V 

F IN5 R OUT3 Z GND 

H IN6 S OUT4 

J IN7 ue OUTS5 
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4.2 Analog-to-Digital Conversion 

An Intel 1H4108, eight-channel multiplexer is used to 
pass the different signals to a Teledyne Semiconductor 
8700CJ ADC. The 8700 is fairly typical of numerous ADC 
devices currently available. It is an 8-bit CMOS chip housed 
in a 24 pin package, capable of up to 800 conversions per 
second. Cost was about $15. Only a negative reference 
voltage, dual power supply, and a few passive support 
elements are required for operation. Any desired positive 
input range can be selected (nominally 0 to 5 V here). The 
circuit, which closely follows the recommendations of the 
manufacturer, is shown in Figure 4.3. The positive power 
Supply iS provided by a pencerd 5S volt regulator (7805). 
Negative 5 volts is obtained from this with a unit-gain 
inverting operational amplifier (741) circuit. An LM336 
voltage reference chip (+2.499 V for 0°C to 70°C) and an 
inverting operational amplifier provide an adjustable 
negative voltage reference. 

There are several advantages to chosing an eight bit 
device. They are much cheaper than 10- or 12-bit units of 
the same speed and they interface easily with an 8-bit 
microcomputer. Byte-size words can be manipulated much more 
etficientlyveand stored, in bess space than longer ones. The 
drawback is resolution. The full-scale resolution of an 
8-bit ADC is 0.4%. While this sounds respectable, if the 
signal fluctuations only cover one tenth of the available 


range then the resolution is at best 4%. For the measurement 
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Figure 4.3. 
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of atmospheric turbulence this points out the importance of 
offsetting the large mean values and amplifying the 
fluctuations so that they ultilize the available digitizing 
range as fully as possible. If this is done properly there 
would be little benefit, given sensor limitations, in using 
more bits. 

The I/O port of the VIP-711 makes interfacing the ADC 
and MUX very simple. The complete circuit is shown in 
Figure 4.3. The output of the ADC connects. directly into the 
parallel input port (IN7-INO). The 8-bit output port 
provides the control signals required. Its three least 
Significant bits (OUT2, OUT1, OUTO) select which multiplexer 
channel is to be sampled. Binary settings of 000 through 101 
are used to designate channels one to six respectively. The 
OUT5 bit instructs the ADC to commence a conversion when set 
high. The MSB of the output port byte, OUT7, controls the 
tape recorder via a low-power relay. When this bit is low 
(the usual state) the cassette recorder is on. When the OUT7 
is high, the cassette is off and the relay draws about 
100 mW. 

The two external flags available on the I/O port are 
also used. EF3 is connected to the external clock circuit in 
Figure 4.3 which is used to set the sampling frequency. EF4 
is connected to the BUSY pin of the ADC thus enabling the 


computer to keep track of the status of the conversion 


process. 


to eonsdsegmt» silt ‘quo! o'sshbed eat ey ‘pa 

arts phiyi: tiga’ bts. seu ley Saaelaioi h +88. 
onisiziaih sidelterd a3 ost t bade ‘yetta dedys an iy 
axa c(lisqoxg, €00b ef ging 3h Lsisteqee te dxion 


pHiey «a ,ehoftaetiarl _gomase: are Sean aEAREE 


= ( 


MA sds pn sosisssre eae Nie ay and. 30 eam) ok <2 
ni Gwode ab. tiuds et vai oar. Meee as Ps 

aft, Ofni, nlgoetie) asosinos ae: ty ae Lucente Ee 
tnod) Seqwuo. trd-8 sit itgae: “ea Jeoke sgt 
1955] #80nd esl |. beatieet 2 feonigha \outniss sits 
eiete lid dyidw Josten (orhd seg areay: akin 
FOr menor s7006. (to epnssisee wigs tines ae ad ag i : _ 
od? wvlevirasaees' vba ot and efsaiene satenpiaeb ot & 
tez “edi. fn betevncs 5 Sonbminob oF OK ‘end gsauivent | 


sii? atodiyps:< FIMOM Baya 3984 pee. “one te. eee] sie 


Dae i 


vol ay aid"atda ngaw . sles saws & iy, caa wt 
PIGS svt nani, 1f yebtosws 21iteaag PR (aap Sonne 
suds ewe es Sd Ag 29 3 mi adconnet ne, tei | 

eee fica: ree i sees Biers Paniorts owl oat ae ni 
itt riyssbe wens famteds ais G4. tassphs Bit! Stn bom. 


$43 seesntsies’ ae at 3 o2-4 | 


43 


A flowchart of the digitizing process is given in 
Figure 4.4, The collection process begins when EF3 detects 
the rising edge of a clock pulse. This initiates a 
conversion scan intended to simulate simultaneous sampling. 
First, an #80 (binary 1000 0000) is placed on the output 
port which selects the first channel but does not initiate a 
conversion. After a short delay (18 us) to allow the 
multiplexer eve to settle down, an #A0 (binary 1010 0000) 
starts the conversion. Both the #80 and #A0 keep the 
cassette recorder power turned off. Once the conversion is 
Started, the VIP-711 continuously polls EF4 which is 
connected to the ADC Busy pin. When the conversion is 
complete EF4 is set high by the ADC and the computer 
immediately latches up the byte on the input port and stores 
it in the next available memory location. The next 
multiplexer code, normally #81, is then sent out followed, 
as before, by #A1. This cycle continues until all active 
channels have been sampled at which point the computer 
resumes polling EF3. The time per conversion is 
approximately 1.25 ms, hence a six channel scan takes about 
10¢ms,. This sean order of magnitude faster than=the typical 
sampling rate of 10 Hz thus providing an adequate 
approximation to simultaneous sampling. 

The data collection proceeds as above until the RAM 
buffer is filled. The first two pages of RAM (#0000 to 
#01FF) are used to hold the data collection software. The 80 


bytes at the top of RAM are used by the operating system, 
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START e 
INPUT #CHAN ss 
ADC 


BUSY? 
NO 
INITIALIZE ADC 
MUX COMMANDS INPUT BYTE FROM ADC 
AND STORE IN BUFFER 


NO 


YES 
TURN ON 
RECORDER 
READ DATA ‘TO TAPE 
AND RESET BUFFER 
TURN OFF 
RECORDER 


Figure 4.4. Flowchart of analog to digital conversion. 
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hence the last page of -RAM is not included as part of the 
buffer. This leaves a total of 13 pages (3328 data values) 
for the buffer which, of course, must be divided among the 
number of active channels. This can be specified in the 
program, as in the software in the Appendix, entered as 


input via the keypad for greater flexibility. 


4.3 Writing to Tape 

When the available semiconducter RAM is filled, the 
collection process is halted and data is transfered to tape. 
Sampling four sensors at 10 Hz will fill up about eight 
pages of RAM per minute so this dumping must occur quite 
frequently. 

By uSing the subroutines contained in the RCA ROM based 
Operating System the data collection process was implemented 
with only moderate software development. In order to do 
this, however, it was first necessary to disassemble and 
decipher the VIP-711 Operating System. A memory map of this 
software is given in Table 4.2. A documented listing of the 
VIP cassette write software is given in the Appendix. 
Basically; the software, via°Q;” generates one cycle of a 
squaréiwave at“2kHz to représent aelow-bit fandfonetcycletat 
G8 ukhztfor avhigh=bit /'Thedinteriace «hardwareeroundspenese 
into crude sines waves for the audio tape recorder. 

In the data collection software, R3 is initialized to 
the address (#8091) thus allowing the RCA write-to-tape 


routine to be called by making R3 the program counter. In 
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| TABLE 4.2. RCA OPERATING SYSTEM MEMORY MAP 
SS ———— 


Address Function Starting at given address 

8000 Locate top of RAM (from #01FF) 

8022 Run program or enter operating system monitor 
8027 Operating system intializations 

806E Keypad select operating system function 
8091 Tape write routine 

80C2 Tape read routine 

80EF Memory read routine 

80F3 Memory write routine 

8110 Digit display for video output 

8143 Video interrupt routine 

816F Subroutine to write a bit to tape 

8182 Subroutine to read a bit from tape 

8194 Hex keypad input subroutine 


addition, RC is initialized to #816F, the starting address 
of the the RCA subroutine to read single bits to tape. This 
Subroutine terminates with an unconditional branch to 
location #816E where P is reset to R3 thus returning control 
to the exit point in the tape-write rountine. Doing this 
from #816F leaves RC with the starting address of the 
subroutine (816E+1) for future calls. This is the most 
common method of calling COSMAC subroutines. The tape-write 
routine is not, however, such a Subroutine as it terminates 
in a closed loop; thus, an eScape route is needed in order 
to use this software. In addition, values normally supplied 
manually via the hexpad must be generated in software. 
Borbunacely, maustupriormto (falling into this Toop, a video 
display subroutine is called by making R4 the program 


counter. Since a video display is not needed during the data 
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collection process, this subroutine call can be used as way 
out of the operating system. By making R4 the main program 
counter of the data collection program the write routine 
will return control almost like an ordinary subroutine. The 
only difference is that RC is left pointing to the last line 
Of khe-routine,’ ratherithan.the filrst,%after “each cal? v Thus 
R3 must be reinitialized before each call, a small price to 
pay in return for avoiding the need to reproduce the tape 


interface software in RAM. 


4.4 Field Operation 

The VIP-711 comes with a light, plastic case which is 
not suitable for outdoor operation. Figure 4.5 shows the 
plexiglass enclosure constructed to house the VIP-711, tape 
recorder, voltage regulators, RF modulator and relay 
Gimcuit,): In the field, this wnt isits om top of the APUs. the 
two units are connected by a single ribbon cable. An 
external 12 V lantern battery must also be connected to the 
microcomputer unit (MCU). During initial experiments, the 
MCU was allowed to draw power from the APU. This was found 
to be unadvisable, however, since the power consumption of 
the computer (about 3 W) is much greater than that of the 
APU. ‘Ini facts: to extend battery) life, thet power to thes MCU 
should not be connected until the analog gain and offset 
adjustments have been made. 

When ready to begin data collection, the data 


collection software must be loaded from cassette tape into 
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the first two pages of RAM. This is done by placing the 
cassette into the tape recorder and moving to the start of 
the audio code. The VIP-711 operating system is then 
selected by setting the microcomputer to RUN while holding 
the C key depressed. To load a two page program into memory 
Starting at location #0000, the sequence of keys 0000B2 are 
pressed and the recorder is turned to PLAY. The tape light 
goes out when the program is loaded and the software 
cassette should be removed and replaced with a blank 
cassette for data collection. To run the program the 
microcomputer is RESET then switched to RUN. 

Numerous versions of data collection software were 
experimented with; the program given in the Appendix was 
found to be the most useful. Upon running the program, 
execution pauses to wait for two input parameters. The first 
is the number of channels to be used, taken sequentially 
from T, W, U, V, X, and Y. While this may seem somewhat 
restrictive, it ensures no confusion about the contents of a 
tape during playback. The second input parameter is the 
number of pages per channel to be allotted for buffering. 
The software checks both inputs for illegal values. 

Once the two parameters are Supplied from the hexpad, 
the: program»waits for the APU clock to’ bejswitched:to;the 
desired sampling frequency. At this point sampling commences 
and the tape recorder should be turned to RECORD (power is 
off). The microcomputer fills the buffer, turns on the 


recorder, writes to tape, turns the recorder off and repeats 
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this cycle until the computer is RESET. Thirty minutes of 
tape will hold about 200,000 bytes, which typically 
FPepresents® 2-066" 34 hours. of data’ 

The only difficulties encountered during the testing of 
the prototype system were caused by low batteries. Symptoms 
of this occurence included software crashes and failure of 
the relay to operate correctly. New 12 V lantern batteries 
lasted about 15 hours. When using old batteries, the system 
must be monitored closely. An audio signal which is 
generated by the VIP-711 during tape transfers proved to be 


a useful indicator of proper operation. 


4.5 Continuous Data Collection 

The major shortcoming of the prototype data collection 
scheme, gaps in the data which occur during transfers to 
tape, could be eliminated with improved software. The only 
difficulty involves the timing done by the microprocessor as 
it writes individual bits to tape; the microprocessor can 
not be disturbed during this process which makes it 
@ifficult to use the COSMAC interrupt capabilities. However, 
Since it takes about 1 ms to write @ single.bit to tape, it 
would be possible to poll the flag connected to the clock a 
thousand times a second to determine if an A-to-D conversion 
must be started. Since the clock pulse stays high for more 
than 1 ms, is should be possible to implement continuous 
data collection with the prototype system by simply modifing 


the collection program. It would become necessary, however, 
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to include the tape write routine in RAM although the RCA 
subroutine which writes individual bits to tape could 
probably be used directly from ROM. 

The envisioned software would use two buffers, say, one 
page per channel. Data would be stored in buffer A until it 
is filled then into buffer B. While B fills, buffer A would 
be written to tape. Clearly, if no data are to be lost, the 
Sampling rate iS constrained by the requirement that one 
buffer can be written to tape more quickly than the other 
can be filled. Using one page per channel and sampling at 
10 Hz gives 25.6 seconds. For a dedicated tape-write this is 
enough time to save eight pages (2 kilobytes) of data. 
Assuming the modified write-to-tape routine is not 
Significantly slower than the present implementation, it 
should be possible to collect six or seven channels of 


continuous data iat. 10 Hz. 
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Chapter 5 


TRANSFERRING DATA TO A HOST COMPUTER 


5.1 A Serial interface 

The data stored on cassette tape must be transferred to 
a host computer before any serious processing can begin. The 
host computer available at the University of Alberta is an 
AMDAHL 5800/5860 operating under the Michigan Terminal System 
(MTS). The University's Computing Services supports a large 
number of remote terminals connected to the computer by 
modems (Gandalf Data Sets) capable of serial, full duplex 
data transmission at speeds ranging from 110 to 9600 baud. A 
Decwriter terminal, normally operating at 600 baud, was 
available for use in this project. The connection between 
the terminal and its modem is via the widely used RS-232C 
Standard (see Figure 5.1) which defines interconnect 
Signals, voltages, connector type and pin assignments. 
Information sent to the AMDAHL must be ASCII (American 
Standard Code for Information Interchange), with odd parity. 
Thus, in corder. for the VIP-7 11> to use the™serial 
communication link, data must be retrieved from cassette 
tape, translated from bance to ASCII code and transmitted, 
Onemputeat actime, in RS~232C compatible Sugqnals. Almost; ell 
of this could be done in software. The only hardware 
additions concerned the RS-232C hook up. 

The RS-232C interface requires that a false bit be in 
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Subminiature D Connector 


1.Protective Ground 
2.Transmitted Data 
3.Received Data 
4,.Request to Send 
5.Clear to Send 
6.Data Set Ready 
7.Signal Ground 
8.Data Carrier Detect 
9...14.Not Used 


15.Transmitted Bit Clock Internal 
16.Not Used 

17.Received Bit Clock 

18,19.Not Used 

20.Data Terminal Ready 

21.Not Used 

22.Ring Indicator 

23.Data Signal Rate Selector 
24.Transmitted Bit Clock 

25.Not Used. 


Figure 5.1. The RS-232C standard. The pin definitions above 
refer to the terminal side of the connector. In the present 
application, only pins 2, ed,eanawy ane of interests Athigh 
Dire must pe 1n bhegrange -S tom-i5.V > whi lekawlow, bits 
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-5 V to -15 V. The buffered Q line on the VIP-711 I/O port, 
used to produce the serial output, generates 0 or 5 V_ 
Signals. While special integrated circuits are available to 
transform such signals to RS-232C levels, the conversion can 
be easily done with more common components. The noninverting 
Operational amplifier configuration shown in Figure 5.2 was 
Msed here. The circuit (transforms 0: V to =10 V and 5 V’ to 

10 V which seems to be the opposite of what is required. 
This is done because the AMDAHL requires a high signal when 
no data is being sent but Q is normally low. Using the 
Circuit shown, Q must be turned on to represent a false bit 
and off for a true bit. This is easily done in software. 

The transformed serial output must be connected to the 
receive line on the modem. However, rather than simply 
replacing the normal terminal with the VIP-711, the 
connections shown in Figure 5.2 leave the usual terminal 
host connections intact while allowing the VIP-711 to tap 
Picomonednop outsof tCheycommunmcationy link. Thisvellows, tor 
example, the standard terminal to be used to log onto the 
system before the VIP-711 is engaged to send data. At the 
end of a data dump, control can be easily passed back to the 
terminal ‘to allow manipulation of the data file without 
@isconnecting the VIP-711. The interface switches on the 
front panel of the microcomputer enclosure (see Figure 4.9) 
allow transmission of information from the microcomputer to 


the terminal and the host as well as the normal terminal to 
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Figure 5.2. Hardware interface between the VIP-711 and a 
host computer and terminal. The serial interface was oper- 
ated at 600 baud to allow simultaneous transmission to both 
the terminal and the host computer. 
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5.2 Reading from Tape 

Reading data from the cassette under program control is 
Similar to the writing of data previously discussed. The VIP 
hardware roughly converts the audio signal back into a 
Square wave which is monitored by EF2. A bit begins when EF2 
goes high and the VIP software determines whether the bit is 
true or false by timing from the rising edge to the falling 
edge. R3 serves as the program counter for the read routine 
which begins at #80C2. This time RC must be intialized to 
#8182, the subroutine which reads in individual bits. Again 
R4 is used as the main program counter so that it is 
possible to escape from the operating system. The total 
number of pages to be read from RAM must be placed in RE.0O. 
In the software given in the Appendix, this value is 
calculated from the number of channels and pages per channel 
which are supplied from the hexpad along with the number of 
channel scans to be printed per line. The output buffer 


begins at location #0200. 


5.3 Binary to ASCII Conversion 

After filling the RAM buffer the VIP-711 must begin 
sequentially sending ASCII data to the AMDAHL. ASCII is 
designed ho handle alphanumeric information and, as a 
result, is rather inefficient for transmitting pure numeric 
data as must be done here. For example, the decimal number 
234 is represented by the single binary byte 1110 1010 


(#PA). To transmit this value in ASCII requires a byte: for 
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each decimal character plus a byte for the ASCII space 
character. 

The conversion from binary to ASCII is aided by the way 
the ASCII code is set up. Adding 30 to a number between 0 
and 9 gives the resulting ASCII hex code. Therefore, the 
problem is to break the binary representation down into its 
three decimal digits. The procedure, depicted in the flow 
diagram in Figure 5.3, begins by subtracting decimal 100 
(#64) from the data byte. If the remainder is less than 
zero, the first digit is zere. Ifjnot 4an additional 
Subtraction is performedyto.find out if this digit is a sone 
Or a two. The remainder of the last successful subtraction 
is then subjected to a generalized version of this procedure 
where 10 (#0A) is successively subtracted until the tens 
Gigit is found. The ones digit is automatically produced as 


the remainder of the the tens loop. 


5.4 Parallel to Serial Conversion 

As each ASCII character is generated it must be 
Serlally transmitteanto the Moses In addation to thesseven 
ASCII bits*associated-with-each digit, a-start—bit)..odd 
parity bit, and stop bit must*bejtransmitted einexpensive 
monolithic integrated circuits called UART's (Universal 
Asynchronous Receiver Transmitters) are available to provide 
this formatting. Here, however, it was implemented in 


eortwarevat no. extra cost. 
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Figure 5.3. Flowchart for Binary to ASCII conversion. The 
Flowchart starts with the binary data byte in the data (D) 
register. The hundreds (HUNS), tens, and ones digits are 
sequentially determined and transmitted. 
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A true bit is generated by setting Q low (false bit 
requires Q high) and falling into a software delay loop. For 
600 baud the necessary delay is about 1.6 ms.'° In the 
actual software, the time of the delay involves some 
overhead due to additional instructions executed before and 
after entering the loop. Since the number of these 
additional instructions varies depending on where the delay 
Subroutine is called from, the number of iterations. is sent 
as an input parameter to the delay subroutine. Values were 
determined with the aid of an oscilloscope, however, the 
setting was not particularly critical. Fine tuning would be 
more important for higher baud rates. It seems likely that 
1200 or 2400 baud transfers could be nprencntee using the 
delay loop approach. 

The serial transmission of a digit, flowcharted in 
figure 5.4, begins by sending a low start bit to the AMDAHL. 
Then, with the data byte in the D register, each bit of the 
data byte is sequentially shifted into DF. By testing to see 
if DF is a one or zero, Q can be set to the corresponding 
value. One of the scratchpad registers is used to generate 
fhe parity bitewtinvralty zero, each time athagn™ bic Ts 
encountered this register is incremented by one. When all 
seven bits have been transmitted, the LSB of this register 


‘©The simplest COSMAC delay loop has the following form; 
loop DEC Rn; Decrement Rn 
GLO Rn; Place Rn.O in the D register 
BNZ: /OOpsgr colo) J Oop.1t. Rn..0= 0 
where the intial value of Rn determines the delay time. Each 
of the instructions requires about 9 us to execute giving a 
27° us delay per loop. This requires, 60 (#3C) repetitions for 


1.6 ms. 
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START © 
INITIALIZE: 
: DELAY ROUT 


PARITY<--0 


COUNTER<--COUNTER- 1 
SET Q 
| NO IS 
DELAY ROUTINE ' COUNTER 
=0? 


YES 


SHiWl Deat© LAE 
RIGHD BYe dB TT 


DF<--D(MSB) 
D(C LSE) <--0 


YES 
RESET Q SET Q 
{ DELAY nourine | 
RESET Q 


RESET Q 
PARITY<--PARITY+1 
i DELAY ROUTINE : 


Figure 5.4. Flowchart for parallel to serial conversion. 
Pereaghdgh bithdDEFi)entens@necessary to” reset’ O. For a’ low 
bite \(DE=0) 00) 1sset ion. 
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will contain the required odd parity bit. This is 


transmitted and followed by a high stop bit. 


5.5 Operation 

With the microcomputer connected into the terminal 
modem link, and interface switches set at T= TRM and 
R= HOST, the user logs onto the host system in the usual 
manner. For the AMDAHL/MTS system used with the prototype, 
it 1S important to issue a %PRENT=ON command before- 
attempting any data transfers. This command activates input 
line buffering by the front end processor of the host 
system. If this feature is not used, occasional characters 
will be lost when the computer, due to its time sharing 
Structure, 1S unable to respond immediately. 

The serial interface program is loaded from tape (same 
procedure as used to load the data collection software in 
Section 4.5) and the data tape is placed in the recorder. 
The microcomputer is Switched to RUN. In the prototype 
software, the number of channels, number of pages per 
channel, and the number of channel scans to be stored per 
line of the host computer file must be entered via the 
hexoad. vAt thisspoint the Poswitohiis set {te COS (COSMAG), 
the cassette is turned to PLAY, and the last key pressed is 
pressed again to initiate the transfer process. 

At 600 baud, the transfer of data is rather slow. The 
conversion from binary to ASCII results in four times more 


characters to be transmitted than were originally saved. 
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However, it is convenient to use 600 baud, since it allows 
easy switching between the VIP-711 and the terminal, as well 
as a Simultaneous hard copy output. Parity errors are 
occasionally detected when reading data from tape to RAM. 
When this occurs, the VIP-711 link with the computer is 
automatically broken. It is then necessary to rewind the 
tape to the start of the bad block and reread it. In many 
cases the parity error does not reoccur. If the error does 
persist, it can be ignored (likely just a single bit) or the 


entire block can be skipped. 
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Chapter 6 


DATA ANALYSIS 


6.1 Experiments 

Considerable turbulence data were collected during a 
series of initial experiments intended primarily to test the 
performance of the various components of the prototype data 
acquisition system. High-level processing was performed on 
some of these data in order to compare computed values with 
those found by others over flat, uniform terrain thus 
providing an additional, although limited, check on the 
operation of the system. 

Most of the experiments were conducted at the 
University of Alberta Research Farm located about 10 km 
south of Edmonton. A permanent, open constuction, 10 m tower 
is located at this site which is a modest approximation to 
horizontal homogeneity. Experiments considered here were 
restricted to days when the wind was predominantely from the 
the west, as this provided the most uniform, unobstructed 
approach to the tower. Even in this direction, however, the 
land sloped downwards at a rate of about one in thirty, and 
hadital Eetchwof tonly about3s00rms Ground cover arouna™the 
tower was short grass. About 50 m to the west this changed 
to a barley crop which was 30 cm high in early experiments, 
but cut to short stubble in later ones (those in October). 

Typically, turbulence measurements were made at about 


5 m. Attempts were also made to obtain supporting wind and 
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temperature measurements; however, the difficulty of 
obtaining reliable profile measurements should not be under- 
estimated. In most cases, the available equipment for these 
measurements, Rimco impulse cup anemometers and fine-wire 
chromel-constantan thermocouples, were inadequate for 
resolving the small gradients encountered. Consequently, 
profile estimates of stability and surface layer parameters 
are quite crude. The poor profile data, site imperfections, 
and the low frequency cut-off imposed by the prototype data 
acquisition software make it unproductive to compare profile 
and covariance flux estimates. Here, therefore, the focus is 
on the spectral analysis of turbulence which the collection 


of raw data makes possible. 


6.2 Sonic and Platinum wire Power Spectra 

The collection of data in segments containing multiples 
of 256 points is well suited to spectral analysis via the 
Fast Fourier Transform (FFT) method which requires time 
series of length 2” (m an integer). Typically, transforms 
were performed on segments containing L= 256, 512, or 1024 
raw data points. This length determines the lowest frequency 
which can be reliably resolved (1/LAt where At is the 
Sampling period) by the subsequent analysis. Each segment 
was detrended with a least-square-fit and its ends (10%) 
were tapered with cosine bells prior to transforming as 
described by Bingham et al. (1967). The complex Fourier 


coefficients were calculated with a standard FFT program 
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then squared, summed, and scaled to provide modified 
peridograms.'' Several consecutive peridograms were averaged 
to provide a relatively stable estimate of the true 
spectrum. 

Following the usual convention in micrometeorology, 
loginS(n)} is plotted against log £. Such a graph, unlike 
one of log{S(n)} vs log f, maintains an area-variance 
equivalence. Figure 6.1 shows a typical vertical velocity 
Spectrum plotted in both ways. The spectrum is the result of 
averaging 20 consecutive peridograms covering a period of 
about 15 minutes. While the large scatter implies the need 
for longer averaging periods, the strong influence of 
diurnal variations restricts such an approach; rarely can 
stationarity be assumed for more than 30 minutes in the 
atmosphere. Following current practise, a smoother picture 
is obtained by block averaging the spectral estimates over 
logarithmic intervals to give about 20 equally spaced values 
on a logarithmic frequency scale. On such a graph, the high 
frequency values involve averaging over many more estimates 
than do lower values and, therefore, are more stable. In 
Darticular™ not too much contwdence can ber placed iingthe 
first few low-frequency estimates. 

Figure 6.2 shows the block-averaged logarithmic spectra 
of sonic vertical wind, nSw(n), and platinum wire 
‘ithe standard FFT routine distributes the variance over 
L=2™ estimates, of which only the first L/2 are unique. 
Therefore, the spectra are scaled by multipling by 2 and 
dividing by An such that o7= Z4S,An in keeping with the 


definition for power spectra used in micrometeorology, 
o*= foS(n)dn. 
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Figure 6.7. A typical vertical wind spectrum obtained® from 
15 minutes of 10 Hz sonic anemometer data. Twenty segments, 
each containing 256 values were Fast Fourier transformed to 
give periodograms with 128 estimates, each which were 
averaged together. The lower plot (diamonds) is log{Sw(n) } 
versus log f and therefore is expected to follow a -5/3 
Slope in the inertial subrange. The upper plot (squares) is 
the logarithmic spectra, log{nSw(n)} versus log f and should 
Ghus, i£0)4 ow aya2y¥eis lope - 
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Figure 6.2. Block-averaged sonic w-spectrum (a) and 
platinum-wire @-spectrum (b) from 30 minutes of 10 Hz data 
obtained on October 10, 1982. Solid lines are from Kansas, 
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temperature, nSg(n), from 30 minutes of turbulence data 
obtained at the Ellerslie site between 1330 and 1400 local 
time during an experiment on October 10, 1982. The spectra 
are the average of 20 modified peridograms, each obtained by 
transforming time series consisting of 512 points of 10 Hz 
data. The wind speed at 5 m during the measurements was 
about 5 m/s, thus making the effect of vertical line 
averaging negligible. The data were low-pass filtered at 
5 Hz; thus, as discussed in Chapter 3, no spectral 
corrections were applied. For these, and all following 
measurements, the sonic anemometer was only calibrated 
relative to a Gill propeller and may, therefore, be in error 
by +10%. A more rigorous wind-tunnel calibration is required 
before the sonic anemometer is used for more quantitative 
research. Tilt and flow distortions caused by the sensors 
introduce additional uncertainties which are difficult to 
quantify. 

The solid lines superimposed on the spectra in Figure 
6.2 are from the extensive measurements by Kaimal et al. 
(1972) made at a near ideal site in Kansas. Kaimal et al. 
normalized their spectra such that all stabilities would 
eOnverge tofa single"eurvesin*thevinertralesubrange. Here, 
the lines from Kaimal et al. have simply been fitted by eye 
in the well-defined inertial subrange. Shown in the figure 
are the Kansas curves for $= -0.1 and -0.3. Profile data 
obtained during the experiment indicated $= Ri= -0.10, hence 


there is fairly good agreement in the shape of the Ellerslie 
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temperature and vertical wind spectra with those from 
Kansas. Although not too much faith can be placed in the 
first couple of low-frequency estimates, it might be noted 
that a departure of the spectra from the reference curves 
seems to occur around f= 0.1. Given the wind speed of 5 m/s, 
this corresponds to wavelengths of around 50 m, the distance 
to the upwind change in surface roughness. 

The spectra of vertical wind and temperature, as noted 
in Chapter 1, can be used to estimate dissipation rates. In 
terms of the logarithmic spectra, nSw(n), and reduced 


wavenumber f= nz/U, (1.4) may be rewritten as 
nSwhf)= ¢4/3ijeles)G.7(27£)74°. TS 


Using a= 0.55. from Dyer and Hicks (1982), and choosing 
f= 1.5, which is well into the subrange yet still free of 


aliasing, gives 


NSwi ld) =n Ones Cegy oy (6a 29 


Oteat io MM, 


e=uL2n08onSwhive) ds. - b6.34 


From figure 6.2(a), nSw(1.5)=0.0080 m?/s* which gives 


e=0.0022 m*/s?. For near neutral conditions 
=e gy Kae (6:..4)) 


which implies a friction velocity of about 16 cm/s. A 
profile estimate of uy = 20 cm/s was obtained on this 


occasion, indicating that the dissipation estimate is quite 
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reasonable. A covariance estimate can not be given, since 
horizontal propeller data were not recorded during this 
experiment. The value of e, unlike a covariance estimate, 
should not be affected by the low-frequency cut-off imposed 
by processing short segments of data. Implicit, of course, 
is the assumption that stationarity existed for a sufficient 
time prior to the measurements, to allow the inertial 
Subrange to adjust to the present conditions, which was 
likely the case. 

Using B= 0.68, from Dyer ‘and “Hicks” (1982) P"and@f="1.5, 


the dissipation of temperature variance, N, (1.5) becomes 


n$4(1.5)= Owt52 Ne~1/%z2/3, (6.5) 


On at z=5 mM, 


N= 2,25 72 eanSa i>). 6516;) 


From Figure 6.2(b), nSg(1.5)= 0.0038°K*. Using the previous 
value ‘of e,/ gives N= '0.0011°K*/s. For near neutral 


conditions, 


= Ha. CoG auekal, C670 


Za 
2 


thus the subrange value of N corresponds to a heat flux of 
around 25 W/m’. This agrees fairly well with a concurrent 
covariance estimate of 35 W/m* obtained by averaging ten 
blocks, each containing 1536 data values (the standard 
deviation of this estimate was 12 W/m?). Again, therefore, a 


reasonable dissipation rate estimate is obtained from the 


inertial subrange. 
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For the October 10 data, stationarity over the 30 
minute period considered was probably a good assumption. 
There are, however, situations when this would not be the 
case. On ,September--16,,.1982,..the prototype data acquisition 
System was used to collect 15 minute sequences of 10 Hz data 
at 45 minute intervals from 1700 to 2100 local time. During 
this evening transitional period, the nature of turbulence 
undergoes large variations. Figures 6.3 and 6.4 show the 
vertical wind and the temperature spectra obtained at 
P780-1715¢" 1830-1845, and 2000-2015 which™correspond to 
unstable ($= -0.1), near neutral, and stable conditions 
($= +0.1). The wind remained fairly steady during this 
period and a single mean value has been used to compute the 
reduced wavenumber. For the vertical wind spectra there is a 
slight indictation of a relative shift of energy to higher 
frequencies with increasing stability, Superimposed on the 
overall drop in the intensity of turbulence. 

The changes in the magnitude of the temperature spectra 
are more dramatic than those in vertical wind. The intensity 
of the temperature fluctuations drops to a minimum at 
neutral then increases rapidly as the stable stratification 
sets in. The breakdown of inertial subrange behaviour in the 
neutral spectrum, as well as the high frequency endsof) the 
unstable spectrum, is probable due to noise problems which 


are discussed in Section 6.4. 


2 


OE ade 1370) igiininots 
TOLER ES Bodh! a" dds 
attd'.2d 29%) Siaow wit: an 


notele typos esen eye ese za: a aly oe andiniagnt 
Sveh SH Ot “20 zantaupag aH TEM: Sr patios 02 ang. am 
oatid . ani? weeel ow ow Gait: nba? ee aa 
so0sildsy? Po -satchan ala Bobyed! fsholaieneed ¢ 
aid’ kode Ma Bina. fia +9 gt aad lie eief «£ 
$a Ran badds 243 ode suaecigemiee! ae bea baiw he 

oF hndgeenred tothh. or 9gs <i are aber char BEN 

" poor ESnos Sldese Bina fata sis hy Q~ #2) of 


ih 
“ae > ey ae DP 


aint prttye: va dede y late abaieeiy brie iT ste 

ei? ssuqmas (Od) Seay need 2 att 50 fev gon pipnia « bag & 
& si Baents aricece bis Dasha rey wis ih ehiediees 5 
tsfa th oF youens 36 feine ay $e tas, & 46) miseseieeh 5 
due ei tisaze oartmap na leit 
Learaitettad 20 cama Seles lel | 
eronade guiteyggnsa git) 36) abit E vHpuem Gata ot sundiaidiie. 


griensini: ad? ba} ¥ here He. ian 


j org ats aiient s100 


; 


sii). iro Penatiy 


72 


nr 
4) 
S 0.0 
ZF 
= 
.¢p) 
G 
- 
ee 
eH 
O 
fx) 
A, 
Ww 
ps 0.004 
= 
fe) 
ov 
O Unstable ({ =-0.1) 
@ Near neutral 
AuStable? (G00) 
0.0001 
Le ; 
f=nz/U 
Figure 6.3. Sequence of vertical wind spectra derived from 


sonic data obtained during sunset at the Ellerslie site on 
September 16, 1982. Transforms were done on segments of 256 
values. Each spectrum is the average of 20 periodograms 
covering about 15 minutes. Filter corrections have been 
applied and aliasing corrected for by interpolating the 
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Figure 6.4. Sequence of platinum wire temperature spectra 
obtained during sunset at the Ellerslie site on September 
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6.3 Propeller Spectra 

The quality of the previous temperature and vertical 
wind spectra becomes clear upon examining the spectra 
obtained from propeller anemometers. An experiment on 
October 9, under near neutral conditions, included 
measurements of a vertical propeller mounted close to the 
sonic anemometer. Figure 6.5 shows a comparison of the 
vertical wind spectra of these two instruments obtained from 
30 minutes of 10 Hz data. The failing response of the 
propeller is clearly seen above f= 0.3. Figure 6.6 shows the 
phase and coherency spectra between the two sensors which, 
if the response of the sonic is perfect, can be interpreted 
as the transfer and phase functions of the vertical 
propeller. The solid lines superimposed on the data are 
obtained from the horizontal frequency response model given 
im 2.3) pana (2.4) assuming, “arbitarily, ‘that m= 1.5 ms The 
comparison is encouraging, with the exception of the phase 
function at high frequencies. Although the phase estimates 
become: less reliable as) the coherenysdrops off, the return 
to zero behaviour of the computed phase is, in part, 
probably due to the physical Separation Of athe two Sensors. 
The dashed line in Figure 6.6(b) is the modelled phase after 
making a simple correction, based on Taylor s hypothesis, 
for a 20 cm along-wind separation. 

The response of properly aligned horizontal propellers 
is expected to be somewhat better than that of a vertical 


propeller, although no simple method of verifing the 
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Figure 6.5. Comparison of a Sonic w-spectrum and a 


vertical Gill spectrum. Spectra are computed from about 30 
minutes of 10 Hz data obtained on October 8, 1982. 
Transforms were on segments of 512 values. No spectral 


corrections have been applied. 
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transfer function given in Chapter 2 was available. Figure 
6.7 shows the u-spectra for 30 minutes of data concurrent 
with the vertical wind data used in Figure 6.6. The 
longitudinal component was obtained by cosine correcting the 
data from the orthogonal pair of propellers, one of which 
was closely aligned with the mean wind. The-solid lines from 
Kaimal et al. (1972) have been fitted in the region around 
f= 0.2 where the propeller response correction is small, yet 
the gtati seen confidence is still fairly high. The 
placement of the line in Figure 6.7 implies, by comparison 
Wren Kaimal et alt (1972), amvalue of u, of about 20 cm/s, 
which is in good agreement with concurrent profile and 
covariance estimates obtained on this occasion. As expected, 
however, the uncorrected spectrum falls off much too quickly 
at high frequencies. 

When the horizontal propeller spectrum is corrected, 
the high frequency end shoots up sharply in an aliased 
fashion. This apparent aliasing, also seen in the vertical 
propeller spectrum shown in Figure 6.6, was typical of all 
computed propeller spectra. However, the sharp high 
frequency cutoff of these sensors should result in a 
reduction of aliasing. The observed effect, therefore, is 
probably due to noise rather than aliasing. The source of 


this noise is probably the digitizing process. 
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Figure 6.7. Longitudinal wind spectrum computed from 30 
minutes of 10 Hz propeller data. The circles indicate the 
spectrum computed from cosine-corrected data. The squares 
indicate the spectrum after attempting to correct for 
propeller frequency response using L= 1.1 Mm. 
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6.4 Noise due to ADC Resolution 

The validity of using an 8-bit ADC depends upon the 
severity of resolution related problems. The error 
introduced by the digitizing process is bounded between +R/2 
where R is the resolution of the ADC. It seems reasonable to 
assume that, for turbulence data, any value within the 
allowable range is equally likely and the mean will be zero. 
This corresponds to a probability distribution function 
which 1S a constant 1/R between +R/2 and zero everywhere 


else. The variance of such a distribution is 


+R/2 
ont=f (1/R)x?dx= R?/12 (6.8) 
—R/2 


If the actual turbulence data are statistically independent 
of the error introduced by digitizing and this noise is 
white then there is expected to be a contribution Of ER i 2a, 
to each spectral estimate. By reducing the resolution of 
actual data to simulate the effects of an A-to-D convertor, 
this has been shown to be an adequate description of the 
noise problem. As an example, Figure 6.8 shows a 15 minute 
vertical wind spectrum computed from data which had their 
resolution numerically reduced from 0.80 cm/s (8 bits) to 13 
cm/s (4 bits). Despite this large reduction in resolution, 
the effect on the overall variance is only about 10%. The 
solid line on Figure 6.8 shows the ADC error expected on the 


basis of the previous argument and it is in good agreement 


with the computed spectrum. 
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Figure 6.8. Effect of ADC resolution on a typical power 


atmospheric power spectrum. The triangles indicate a 
vertical wind spectrum computed from a set of 8-bit data 
which makes good use of the full range of the ADC. The 
SQuares indicate the spectrum of the same set of data 
computed after reducing the resolution of the data to 4 
bits. The solid line indicates the expected effect of 
reduced resolution assuming it acts to introduce white 


noise. 
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The resolution problem has affected the propeller data 
most severely because relatively low gain settings were 
customarily used in order to accomodate potentially large 
variations which could be caused by wind shifts. Although 
high-pass filtering reduces this problem, it was avoided 
because of the need to reconstruct the absolute signal 
levels in order to perform cosiné corrections. Of course, 
improper use of the gain settings can make resolution a 
problem for any channel. The neutral temperature spectrum in 
Figure 6.4 is probably suffering from resolution noise. 

The ADC noise on any two channels should be 
uncorrelated, thus covariance and cospectra calculations 
Should be fairly insensitive to ADC resolution. Following 
the approach used for power spectra, Figure 6.9 shows the 
cospectra of vertical wind and temperature obtained from 
data with 8 and 4 bits of resolution. The two cospectra 
agree very well. The differences which do occur at high 
frequencies are associated with amplitudes well below the 


noise level. 


6.5 Cospectra 

Cospectra reveal the relative importance of different 
scales of motion to the total flux and are thus of 
considerable interest. Similarity theory predicts an n°’’? 
behaviour in the inertial subrange of both Cwe and Cwu, 


which should show up as a -4/3 slope on the usual 


logarithmic plot. Figure 6.10 Shows nCwe corresponding to 
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Figure 6.9. Effect of ADC resolution on a typical 
atmospheric cospectrum. The line indicates a cospectrum 

nCwo computed from 8-bit data. The circles are the cospectral 
estimates after reducing the resolution to 4 bits. 
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the earlier sonic and platinum wire power spectra for 
October 10 (Figure 6.2). The cospectra exhibit more scatter 
than the corresponding power spectra. This is indicative of 
the need for longer averaging times for cospectra (fluxes) 
relative to power spectra (variances) as discussed by 
Wyngaard (1973). In the inertial Subrange of Figure 6.2, nCwg 
seems to fall off just slightly more rapidly than 
expected. Since the area beneath this curve is proportional 
tomtne-heat flux, the resulting loss of flux at the low 
frequencies is clearly seen. 

It would be nice to be able to show a graph of momentum 
cospectra of the same quality as Figure 6.10. Figure 6.11 
shows nCwy obtained from the October 9 data. It 1S not even 
possible to plot these data in the usual manner, since 
numerous sign changes occur; instead, a semi-logarithmic 
graph of nCwy is given. While the use of horizontal 
propellers is expected to degrade the momentum cospectrum, 
site and actual atmospheric anomolies have probably 
contributed to the poorly behaved result in Figure 6.11. 
Stationarity was expected to be quite good and this is 
confirmed by the agreement between the 15 and 30 minute 
cospectra in Figure 6.11. 

Somewhat more encouraging are the nCwy graphs shown in 
Figure 6.12. These three cospectra were obtained from the 
September 16 data corresponding to the power spectra 
considered in Section 6.2. They have been corrected for 


differences in propeller and sonic response as discussed in 
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Figure 6.10. A cospectrum of vertical wind and temperature 
computed from 30 minutes of sonic and platinum wire data 
corresponding to Fig. 6.2. The solid lines are from Kaimal 
etealen(1972) . 
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Figure 6.11. Momentum cospectra, -nCwy , from 10 Hz data 
obtained on October 9, 1982. The solid line corresponds to 
30 minutes of data while the dashed line was obtained using 
Ohdy using, the fi tet. > minuees No Spectral” corrections 


were applied. 
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10 Hz data obtained on September 12, 


minute segments of 
1982. Corrections for propeller response, 


Chapter 2, were applied. 


as discussed in 
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Chapter 2. At higher frequencies than those shown, the 
spectra change sign and, therefore, were dropped. The 
Scatter of the near neutral nCwy is large, which is not 
Surprising since stationarity is a rather poor assumption 
for this period. Overall, however,these spectra do seem to 
be consistent with a -4/3 slope in the range shown. 
Finally, Figure 6.13 shows the sequence of nCwg for 
September 16. For the unstable Lizus) necessary to. plore 
-nCwe. Again, a lot of scatter can be seen in the near 
neutral cospectrum, but the -4/3 slope of the inertial 
Subrange is easily identified in the stable and unstable 


cases. 
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Figure 6.13. Sequence of heat cospectra, nCyg, from 15 
minute segments of 10 Hz data obtained on September 12, 
1982. In the stable case, -NCwg 1s plotted. 
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Chapter 7 


SUMMARY, CONCLUSIONS, AND RECOMMENDATIONS 


7.1 Summary 

A microcomputer-based, micrometeorological data 
acquisition system has been assembled and successfully used 
to collect and transfer turbulence data to a host computer 
for processing. The system is small, portable (powered by 
common 12 V lantern batteries) and relatively inexpensive. 
Meteorlogical instrumentation consists of a sonic anemometer 
for vertical wind measurements, an orthogonal pair of Gill 
propeller anemometers for horizontal wind measurements, and 
a thin platinum wire (25 um diameter) to measure temperature 
fluctuations. A standard Gill UVW array was used to mount 
all four instruments. Cost of the resulting wind-temperature 
array was about $2500. 

Analog conditioning, muitiplexing;, "and "8-bit analog- 
to-digital circuitry which interface the instruments to an 
REAPVIP=7 11 @singlé-board microcomputer™have beenvassembled 
and described. The microcomputer controls data collection at 
5, 10, 15 or 20 Hz (same for all channels). Simultaneous 
Sampling is approximated by fast sequential digitizing (1.25 
ms per channel). There is enough on-board memory to buffer 
up to 3328 data bytes. When this available Space is filled, 
data collection is stopped and the data are transfered to 
standard audio cassette tape at 100 bytes per second. The 


possiblity of achieving continuous data collection at up to 
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100 bytes per second has been considered and appears 
feasible with improved software. The microcomputer also 
controls the retrieval and transmission of the data to a 
host secomputer via-an RS7232C interface. :iThis-process 
involves reading data from tape, converting them from binary 
Geer oCh scode ,prandythen jtransmittingsathe ASCII bytes in 
Serial form at 600 baud. This interface was implemented 
predominantly in software. The only hardware required was a 
Simple RS-232C driver. Cost of the computer, tape-recorder, 
and electronic components required for this project came to 
around $500. Documented assembly and machine language 
listings of the prototype software are provided. 

A series of experiments to test the prototype system 
were undertaken at the te eee of Alberta Research Farm 
(Ellerslie) located approximately 10 km south of downtown 
Edmonton, Alberta. Statistical and spectral analysis was 
performed on many of the data sets obtained. The results of 
a few selected cases were presented. The -5/3 and -7/3 
behaviours expected in the inertial subrange of atmospheric 
power spectra and cospectra, respectively, were generally 
found to be followed fairly closely. Reasonable estimates of 
dissipation rates were obtained from the subrange of sonic 


vertical wind and platinum-wire temperature spectra. 
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7.2 Conclusions 

This work has demonstrated that it is currently 
possible to construct a microcomputer-based data acquisition 
system suited to the study of atmospheric turbulence for 
around $500 (not including the cost of meteorological 
instrumentation). This is at least an order of magnitude 
less than any comparable commercial systems currently 
available. Hardware assembly is relatively easy and neither 
specialized development tools nor a great deal of technical 
expertise are required. 

The RCA COSMAC VIP-711 training microcomputer was found 
to be well suited to this project, operating reliably during 
two months of intermittent field testing. The commercial 
quality electronic components of this-and other comparable 
microcomputers does, however, restrict exposed operation to 
temperatures above 0°C. 

The analysis of data collected with the prototype 
system has confirmed that the low-cost, CA27 sonic 
anemometer and a very simple platinum wire thermometer can 
be used to study the spectra of turbulence at frequencies up 
to at least 10 Hz and to obtain reasonable subrange 
estimates OfPdissipation rates = Propeller janemometers’, “on 
the otherwhand ,“providerdit tle® intormati onPabove) abouts" “Hz. 
Datayanalysis¢ alsot indicates, that® 8bit analog-to-digital 
conversion can give sufficient resolution for the study of 


turbulence provided that analog amplification and offsets 


are used judiciously. 
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7.3 Improvements 


Being microcomputer-based, the prototype system has a 


great deal of potential for further development. The 


following are a few suggestions for improvements: 


ee 


As discussed in Chapter 4, it should be possible to 
provide continuous data collection at 100 bytes per 
second by modifying just the software. This is the most 
important improvement which needs to be made to the 
CUmnenit. ystems, 

The software could also be improved to permit different 
channels to be sampled at different frequencies. In 
particular, the propeller anemometers could be sampled 
at half the frequency of the sonic anemometer and 
platinum wire thermometer. 

The data collection and RS-232C interface software could 
be transfered to Read-Only-Memory. The VIP-565 EPROM 
programmer and VIP-56C0 EPROM interface board are 
available to facilitate this operation. Placing the 
programs in EPROM would simplify field operation, 
although it should be noted that loading programs from 
tape caused no difficulties during the initival testing. 
Conversion to 12-bit data words would provide ‘sufficient 
resolution to eliminate the need for user adjustable 
gaiineand) offset) control sy) tthuss signif veantly? samplat ing 
fielids operation ofsithe: system. However,icost ofva 
suffaciently fast “2-biteADQis many taumes greater’ than 
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12-bit ADC will be more difficult since the extra 4 bits 
will have to be multiplexed onto the 8-bit input port. 
Using 2 bytes per sample will allow writing only 50 data 
values per second to tape. This should be sufficient, 
however, if continuous data collection is implemented. 
The present 600 baud transfer to the host computer is 
rether slow emoTranster-fates|ofoad200uore 2400ebaud could 
probably be.implemented in software although these 
values would not be particularily useful since they are 
not compatible with the available terminal rates. A 9600 
baud interface which would allow the system to tie in 
with an available CRT terminal would be very convenient 
although this could not be done in software alone; a 
UART chip, baud rate generator, and accompanying 
Circuitry would be required. 

The inclusion of real-time output, notably covariances, 
would be useful in pointing out instrument malfunctions 
Or unusual atmospheric conditions. There is certainly 
sufficent time for the microprocessor to do some 
real-time calculations; in the present software most of 
itse timel ask spentw wad tamge fore ans Apto-Dpconversionnto 
start or finish. Using this spare time, however, would 
require implementing much more elaborate interrupt-based 
software. In addition, since a video monitor is not well 
Suited to field use, a device to provide real-time 
output, preferable a small printer, would also have to 


be purchased and interfaced to the microcomputer. 
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Implementing continous data collection, as discussed in 
Chapter 4, and simultaneous real-time output would 
present a formidable task probably best left to an 


experienced microcomputer system designer. 


7.4 Future Considerations 

While the COSMAC VIP-711 microcomputer seems to have 
been a good choice for this project, RCA has now stopped 
manufacturing it as noted in Chapter 4. In retrospect this 
is not too Surprising. The VIP system was originally 
introduced in 1978; five years is a long time in the micro- 
computer business. There are, however, other training micro- 
computers based on the RCA 1802 processor which could be 
used for this project. The ELF from Neutronics and the 
VENTURE from Quest are two likely candidates. For anyone 
already familar with (or willing to learn) a different 
microprocessor language, there is no need to restrict the 
choice to a 1802-based system. Aside from low-cost, a tape 
Cassette interface running eat 0000 omeamore baud and: a 
sufficient number of decoded input/output lines are the main 
features to look for when selecting a microcomputer. The 
$200 MPF-1 from Multitech, for example, appears to be a 
promising Z8U=based microcomputer, but its cassette 
interface Operates at only T6.5ebytes per esecond which 
makes it unsuitable for the present application. 

There seems to be little hope of seeing reliable 


interfaces for audio cassette recorders at much over 1000 
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baud. Hard disc or bubble memory systems with over a 
megabyte of non-volatile storage have potential, but, at 
present, are still very expensive. Thus, cassette tape is 
likely to remain the cheapest answer to relatively slow, 
sequential data collection for some time to come. 

What about the new generation of 16-bit microcomputers? 
At present, comparable 16-bit single board systems are still 
euate- expensive. |Before™ long, however; ._such systems will 
become widely available and greatly facilitate the use of 
12-bit A-to-D converters in a project such as this. In the 
meanwhile, there is much that can be accomplished with the 


inexpensive 8-bit systems currently available. 
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APPENDIX 


This appendix contains one version of the machine 
language programs used by the prototype data acquisition 
system. The data collection software permits the user to 
input the number of channels, n, to be sampled, and then 
uses the first n channels selected sequentially ErOme Weel 
Uae oo, and Y. lf this order 1s not’ satistactory for a 
particular experiment it can be easily modified by altering 
the initializations in locations #0032 through #0055. 

Both data collection and Host interfacing software use 
routines in the RCA COSMAC VIP-711 operating system. These 
routines have been disassembled and documented and are 
included here for completeness. Except for the RCA 
Subroutines, the programs are loaded from cassette tape into 
RAM locations #0000 to #01FF. Only one program resides in 
RAM at any one time. 

The programs have been developed with the aid of a 
version of the RCA COSMAC single pass cross assembler 
provided by the University of Alberta Computing Services in 
the public file *COSMAC. Certain assembler language commands 
are missing in this limited implementation and machine 
language hexcodes must be used in their place. Most notable 


of these are the REQ (#7A) and SEO (#7B) commands which 


control the Q line. 
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LOC COSMAC CODE 


0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0002 
0003 
0004 
0005 
0007 
0008 
0009 
O000C 
000D 
OOO0E 
0010 
0012 
0013 
0015 
0016 
0017 
0018 
0019 
001A 
001B 
O01C 
001D 
O01E 
OO1F 
0020 
0022 
0023 
0025 
0027 
0028 
002A 
002B 


3011 3 
D7 
AA 
9A 
AF 
94 
AB 
8B 
F4 
AB 
2F 
8F 
3A1B 
8B 
FDOA 
332B 
7B 
3028 
D7 


ENDCH: 
OKCHN : 


PLUS: 


ENDRM: 


oS 


SOURCE LINE 
.»COSMAC DATA COLLECTION SOFTWARE 
. WRITTEN BY PAUL HOPPS (1982) 


- REGISTER DESIGNATIONS: 

eR Pointer to I/O Control Bytes 

--R2.0 Channel Counter 

ae Address of Key-Input 

eeR3 RCA Tape-Write Program Counter 

.-R4 Main Program Counter 

+ oR6 Pointer to next buffer location 
eRe Keypad-Input subroutine P counter 
~-RA.O Number of pages per channel 
»-RA.1 Number of channels 

eRB Total number of pages in buffer 
veRE RCA Write-a-Bit P Counter 

~-RE.O Page counter for RCA routine 

REF Keypad debounce counter 

.RF Remaining buffer space counter 


, #6380 --Turn cassette power off 
GHI 0 »-Initializations, make R4 
PHI 4 ..the main program counter. 
PHP] 

LDI #09 

PLO 4 

SEP 4 

#A1->R7.0 ..Call Key-Input subroutine 
SEP 7 ..to input the number of 
PHI A ..-channels (Store in RA.1). 


SMI #07 -eCheck to make sure the 
BL OKCHN ..number of channels is <7. 
,#7B sePE Sno, PEurn on O*¢eone 


BR ENDCH ..and terminate in a loop. 
SEP 7 .-Else input pages per 
PLO A ..-channel (store in RA.O). 
GHI A .-Get number of channels 
PLO F ..(store in RF.O) 

GHI 4 

PLO B 

GLO B ..-Compute total number of 


ADD .-pages required by 


PLO B »-multiple additions (xX=2) 
DEC F .-Since the COSMAC does 
GLO F »-not have a multiply 


BNZ PLUS ..instruction (store 
GLO weresult. in RB: 0). 
SDI #0A --ls there 

BGE OKRAM ..enough memory? 


tw 


Py alee: SsTEenoF SCurnPon "Oo tone 
BR ENDRM ..and terminate in this loop 
SEP 7 ..else wait for clock to Start. 


F881BC OKRAM:#81->RC.1 
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LOC 

002E 
0031 
0032 
0034 
0035 
0037 
0038 
003A 
003B 
003D 
003E 
0040 
0041 
0043 
0044 
0046 
0047 
0049 
004A 
004C 
004D 
OO04F 
0050 
0052 
0053 
0055 
0056 
0059 
OO5A 
005B 
005C 
005D 
0060 
0061 
0062 
0065 
0066 
0067 
0069 
O06B 
006C 
006D 
O06E 
OO6F 
0071 
0072 
0073 
0074 
0075 
0076 
0078 
0079 
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COSMAC CODE SOURCE LINE 


F8FBA1 #FB->R1.0 

E 1 SEX> 1 

F8A5 LDI #A5 - Multiplexer Initializations;: 
73 STXD 

F885 LDI #85 .-OOFO #80 

rs STXD « OOF 1 -4A0 

F8A4 LDI #A4 ~eQOP Zur #64 

v3 STXD ..O0F Sai FAN 

F884 LDI #84 --OOF4 #82 

Z3 STXD os OOPS) WfA2 

F8A3 LDI #A3 « sU0FG 8453 

nS STXD 00k 7RGAAS 

F883 LDI #83 . OOF By u#84 

73 STXD --O0F9 #A4 

F8A2 LDI #A2 OOF AIT #85 

73 STXD WeUUED esas 

F882 LDIV 482 

7s STXD sin, thelig® control bytes, 
F8A1 LDI #A1 ~--an '8' is used to select 

73 STaBP --a MUX channel and an '‘'A' 
F881 LDI #81 »--1i1S used to initiate an 

oS} STXD acA/D conversion. In this 
F8A0 LDI #A0 seprogreaméchannels 0,1,2,3,4, 
73 STXD sancdwou (Wit ,U,V,%, and Y 
F880 LDI #80 .-respectively) are sampled in 
as STXD sordervupimtootheslast channel. 
F802B6 START: #02->R6.1 ..The buffer 

94 GHI 4 .-Sstarts at location 

A6 PLO 6 ~- #0200 (X=6). 

8A GLO A ..USe RF to keep track of 

BF PHI F .-how much room is remaining 
F8FFAF #FF->RF.O ..in the buffer. 

9A CYCLE:GHI A --Start of a Channel scan. 

A2 PLO 2 

F80FA1 #OF->R1.0 ..Select the first MUX channel 
63 yROS .- before the scan starts. 

21 DEC R1 


3E67 WAITI@BN3CWATIT! ..Wait for a O)to 5 V 
3669 WAIT2:B3 WAIT2 ..transition from the clock. 


63 NEXT: , #63 »--Select the MUX channel. 

C4 ,#C4 ..-Short Pause (NOP) 

63 , #63 .-Commence an A/D conversion. 
E6 SEX 6 

3F6F BUSY:BN4 BUSY ..Wait for end of conversion. 
6B ,#6B --Input the data byte and store 
60 , #60 -e-in the next buffer location 
E 1 SEX 1 

22 DEC 2 .-l1s this the last channel? 
82 GLO 2 .-lf no, then start the 

3A6B BNZ NEXT ..next conversion. 

2F DEC F ..-l1f Yes, check to see if 

OF GHI F Sanne bDutferere.full. 


LOC 
007A 
007¢C 
O07F 
0082 
0083 
0085 
0086 
0089 
008A 
008B 
008C 
008D 
008E 
008F 
0090 
0092 
0093 
0095 
0098 
009B 
009B 
009B 
OOA0 
O0A1 
O0A4 
OO0A7 
00A8 
OOA9 
OOAB 
OOAD 
OOAB 
OOAF 
00BO 
00B2 
00B3 
00B6 
00B7 
00B8 
OOBA 
OOBC 
OOBD 
OOBE 
00C0 


COSMAC CODE 
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SOURCE LINE 


3A60 BNZ CYCLE ..If no, continue the cycle. 
F880B3 #80->R3.1 ..Else get ready for a 
F891A3 #91->R3.0 ..Tape write. 
E4 SEX 4 
6300 , #6300 »-Turn the cassette on. 
EC SEX C 
F802B6 #02->R6.1 ..Tape write 
94 GHI 4 .-initializations. 
A6 PLO 6 
8B GLO B 
AE PLOWE 
D3 SEP 3 .-Call RCA Tape-write. 
6 1 Re oi) JeLurneviTndedc: off, 
E4 SEX 4 
6380 , #6380 -.-Turn the cassette off. 
E 1 SEX (1 
3056 BR FOLART Gta CARTHOVEr . 
000000000000 ,#000000000000 ..Filler 
0000000000 ,#0000000000 
. - KEYPAD) IDNPUT, SUBROUTINE 
D4 RET:SEP 4 ..-Return to Main. 
F800B2 #O00q>R2 WlozkiSetiM(%)s ito 
F8FFA2 #FF->R2.0 ..#00FF (a place to 
E2 SEX 2 ..store value of key). 
92 GHI 2 « Sete-D to; zero 
FCO 1 NEXT2:ADI #01 .-This loop continuously 
FAOF ANI #0F ..-scans the hexkeys until 
52 STR 2 ..one of them is pushed. 
62 #62 
22 DEC 2 
3E6B BN3 NEXT 
7B ,#7B ..-Beep when a key is pressed. 
F820BF #20-> REA 
2F CONT::DEC F .-Debounce delay loop. 
OF GHI F 
3AB6 BNZ CONT 
36BA HERE:B3 HERE --Turn off tone when the key 
7A ,#7A ..is released. 
FO LDX . Recall value of Key to D. 
30A0 BR RET o.- GOCOMNRET:. 
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LOC COSMAC CODE 


0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0000 
0001 
0008 
000B 
000C 
OOOF 
0012 
0013 
0014 
0015 
0017 
0018 
0019 
OO1A 
001B 
001C 
O01E 
OO1F 
0021 
0022 
0023 
0024 
0025 
0026 
0028 
002B 
002E 
0031 


7A 
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SOURCE CODE 


mice Delay loop counter 

ee 2 End-of-line delay P counter 

weR2 ASCII message pointer 

»-R2.0 Character per line counter 

eR RCA Read-Tape-Routine P counter 
.R4 Main program counter 

Asa) Serial Output subroutine P counter 
ed ede Buf ferelogatron fpointer 

Dey: Keyboard Input P counter 

Parity Register 

Number of channels of data 

Total pages in buffer 

ASCII ones digit 

Temporary ASCII character storage 
: Number of characters per line 
ERC RCA Read-a-Bit subroutine 

wCRD Delay routine program counter 

. RE Starting address of buffer 

oeRE Workhorse register 


ww 
> 
—=—s CE © = =1 CE 


rHTA sslnrtialazations: 


SUOBABAZBSBABBBDUORO eilieeR li; R401 )R5. 1,RAL 1, RB. 1,RD. 1 


F80CA4 
D4 
F801B7 
F803A7 
D7 

B8 

D7 
S8AF 
SA 


F827A2 
F881BC 


#0C->R4.0 ..Make R4 the main program 
SEP 4 seCOUNLEN. 
#01->R7.1 
#03->R7 20 
SEP 7 .-Input the number of channels 
PHI 8 .-on the tape (store in R8.1). 
SEP 7 .-Input the number of pages 
RSt 1=2RF.0..<per channeli®von thewtape. 
PAGES:GHI A .-Compute the total number of 
ADD »-pages per block by using 
PHI A »-multiple additions (store 
DEC F ws Desulttain RACIOe 
GLO F 
BNZ PAGES 
SEP 7 .-Input Scans/line to be printed. 
RB 1a>RE-..0 
BYTES: Ghai 7s »-Compute the number of data 
ADD .-values to be stored per 
PHI B .e-line of file in the Host 
DEC F ..computer (store in RB.1). 
GLO F 
BNZ BYTES ..Initializations: 
#AD->R1.0 ..Delay address (#00AD) 
#BD->R5.0 ..P-to-S address (#00BD) 
#27->R2.0 ..Message address (#0127) 
#81->RC.1 ..need by RCA routines 
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LOC 

0034 
0037 
003A 
003C 
003D 
003E 
003F 
0040 
0042 
0044 
0045 
0047 
0048 
0049 
004A 
004C 
004D 
0050 
0051 
0052 
0055 
0058 
OO5B 
005C 
005D 
O05E 
0060 
0063 
0064 
0065 
0066 
0069 
006B 
O06C 
O06E 
0070 
0071 
0073 
0074 
0075 
0077 
0078 
0079 
OO7A 
007D 
OO7E 
007F 
0081 
0082 
0084 
0085 
0086 
0089 
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COSMAC CODE SOURCE LINE 


F8F4AD 
F819AF 
3E3A 
42 

D5 

2F 

8F 
3A3C 
F80D 
D5 
F8B0A 
D5 

D1 

E4 
6300 
EC 
F802B6 
94 

A6 
SAAEBF 
F880B3 
F8C2A3 
D3 

61 

E4 
6380 
F802B6 
94 

A6 

E6 
F8FFAF 
F830 
D5 
OBA2 
F80D 
BS 
F80A 
D5 

D1 
F820 
D5 

72 

AA 
F82FAB 
8A 

AA 
FF64 
1B 
357E 
8B 

D5 
F82FAB 
BA 


HERE: 
MSG13 


LOAD: 


NEXTL : RBs-14>R2. 0 


#F4->RD.0 ..WaitC address (#00F4) 
#19->RF.0 ..25 characters in message 
BN3 HERE ..Wait for repeat of last key 
LDA 2 .. before starting to transmit. 
SEP? 5 -.Ssend the message 

DEC F ~.'$COPY *MSOURCE* -DATA' 

GLO F ..to the Host computer 

BNZ MSG1 

LDI #0D ..Send a carriage return. 

SEP 5 

LDI #0A -.-Send a linefeed. 

SEP 5 

SEP 1 ..-Delay for system response. 
SEX 4 

, #6300 -Turn the recorder on. 

SEX C 

#O025>R61 Vin biaeligations for the 

GHI 4 ~RCA tape read (the buffer 
PLO 6 gascarts at. e200)". 

RA. 1->RE.0,RF. 1 

#80->R3.1 

CAR SU 

SEP 3 ~-Execute RCA tape read. 

,#61 eturn the video off. 

SEX 4 

PHOS p#80 -Turn the recorder off. 
#02->R6.1 

GHI 4 ~-Prepare to send the data. 
PLO 6 

SEX 6 

#FF->RF.O ..RF counts down chars/block. 
LDI #30 -Send a '0' a block sSeparater. 
SEP 5 


~»R2 counts down chars/line. 


LDI #0D ..send a carriage return. 

SEP 5 

LDI #0D ..Send a linefeed. 

SEP 5 

SEP 1 .-Delay for system response. 
NEXTC:LDI #20 

SEP 5 »-send a 'space'. 

pee ..-Get the next data byte. 

PLO A 

#2F->RB.O ..Convert hundreds digit to 

GLO A . ASCII by repeatedly 


LOOPH:PLO A 


..subtracting #64=100 until 


SMI #64 ..the remainder is less than 
INC B --100. Add #30 to the hex 
BDF LOOPH ..result. - 

GLO B 

SEP 5 ..Transmit hundreds digit. 


#2F->RB.O ..Convert tens digit to 
GLO A . ASCII by repeatedly 
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LOC 

OO8A 
008B 
008D 
O08E 
0090 
0091 
0092 
0093 
0095 
0096 
0097 
0098 
0099 
009B 
0O09C 
O09E 
OOA0 
OOA6 
OOAC 
OOAC 
OOAC 
OOAC 
OOAD 
00BO 
00B1 
00B2 
00B4 
OOB6 
OOBC 
OOBC 
OOBC 
OOBC 
OOBD 
OOBE 
OOBF 
00C2 
00C5 
00C6 
00C8 
00C9 
OOCA 
OOCB 
OOCC 
OOCE 
OOCF 
00D1 
00D2 
00D3 
00D5 
00D6 
00D7 
00D8 
OODA 
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COSMAC CODE SOURCE LINE 

AA LOOPT:PLO A . Subtracting #0A=10 

FFOA SMI #0A ..from the remainder to the 

1B INC B .-hundreds loop. 

338A BDF LOOPT 

8B GLO B 

D5 SEPA ~.9ena fens digit. 

8A GLO A ..Ones digit is remainder 

FC30 ADI #30 ,.0f tens 1060p. 

D5 SEP 5 ..Send the ones digit. 

22 DECK 2 

2F DEC F 

OF GHI F .-1f this is the last number 
3242 BZ LOAD ~-in the buffer goto LOAD. 

82 GLO 2 .-Is this the last character 
3A75 BNZ NEXTC...in the line? If no, goto NEXT. 
306C BR NEXTL ..If Yes, goto NEXTL (nextline) 
000000000000 ,#000000000000 ..Filler 


000000000000 ,#000000000000 

~END-OF-LINE DELAY SUBROUTINE 
D4 RETD:SEP 4 .-Return to Main. 
F820B0 DELAY:#20->R0.1 ..Delay for AMDAHL response. 
20 LOOPD:DEC 0 
90 GHI 0 
3AB0 BNZ LOOPD 
30AC BR RETD esGOto RETD. 
00000000000 ,#000000000000 

»-PARALLEL TO SERIAL ROUTINE 
D4 RETS:SEP 4 ..-Return to Main. 
7A PTOS: ,#7A seTRED QOotonesotis 
A8 PLO 8 ~-Binary data value-->R8.0 
F800A7 #00->R7.0 ..Initialize parity reg (odd) 
F807A9 #07.4>R9.0 ..and bit counter (7). 
7B ,#7B ,burn © tones ons 
F83B LDI #3B 
DD SEP D .-Send start bit and then 
88 SHFTR:GLO 8 vestart. Shiftinggoutidata, bits: 
F6 SHR 
A8 PLO 8 .-R8.0 Holds remaining bits. 
33D 1 BDE HIGHK ..Turn Q on for a low bit. 
7B /#7B 
30D3 BR NEXTB ..Get next bit. 
7A HIGH: ,#7A »-PRUrRIO off pionfa kyghebit . 
17 INC 7 .-Increment partiy register. 
F83A NEXTB:LDI #3A 
DD SEP D ..-delay for data bit(600 BAUD). 
29 DEC 9 ..Is this the last data bit? 
89 GLO 9 welteno, get. the next bit 
3AC9 BNZ SHFTR ..else send out the parity bit. 
87 GLO 7 
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LOC 

OODB 
OODC 
O0ODD 
OODF 
OOEO 
O0E2 
00E3 
QOOE5 
00E7 
O0E8 
OOE9 
OOEB 
OOEC 
OOERE 
00F3 
00F3 
00F3 
00F3 
OOF4 
OOFS5 
OOF6 
OOF7 
OOF9 
OOFB 
0102 
0102 
0102 
0102 
0103 
0106 
0109 
010A 
010B 
010D 
010F 
0110 
a 
OnntZ 
0114 
O4sES 
OnaS 
Otho 
011A 
011C 
O11E 
Oil E 
0120 
0122 
0127 
0127 
0127 
O17 
0128 
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SOURCE LINE 


F6 SHR 
A7 BLO 7 
33E2 BDF HPAR 
7B ,#7B 
30E5 BR STOPB 
7A HPAR: ,#7A 
30E5 BR STOPB 
F81B STOPB:LDI #1B 
DD SEP D ~sDetay TOL Darityebpit, 
7A ,#7A 
F8FF LDI #40 s.Send out stop bit. 
DD SEP D 
30BC BR RETS aeReCCUENS 
0000000000 ,#0000000000 
soon RIALsOUTPUT BIT DELAY 
D5 RETC:SEP"5 (4RELUrHNMEO. P-TO-S routane 
AO WAITC:PLO 0 iar tnglembi tade lay Loop. 
20 LOOPC:DEC 0 
80 GLO 0 
3AF5 BNZ LOOPC 
30F3 BR RETC seGOO TEL. 


00000000000000 ,#00000000000000 


D4 RETK : 
FS801B2 INPUT: 
F8FFA2 

E2 

Ce 

FCO! NEXTK 
FAOF 


2F LOOPK: 


367€ PKEY 
7A 

FO 

3002 
0000000000 


24 
43 


meee wee ew ee ae ie ie i = 


» eKEYBOARD INPUT ROUTINE 


SEP 4 ..Return to Main. 

#OK-AR2, 1 seoet M(x) to FOIEF (Ca place 
#FF->R2.0 ..to store the value of key). 
SEX. 2 

GoIenZ 

ADI #01 sefnas Loop Continuously 

ANI #0F ..scans the hex keys until one 


STR 2 ..of them is pressed. 

#62 

DEC. Z 

BN3 NEXTK 

,#7B .»-Beep when a key is pressed. 
#20->RE.< | 

DEC F .»-Debounce delay loop. 

GHI F 

BNZ LOOPK 

B3 PKEY .-Turn off tone when the 
aS ..key is released. 

LDX ssRecCall, value of (key to D 
BR RETK “aula, JO: COmRE TK 
,#0000000000 

. ASCII MESSAGES FOR HOST 

PA2E eS 

HE ST Be G 


a 


mh 
1i¢ Vitis Ter Yatags 


150 .GOle To: DGB e. 


Sti iHos (eevee Os iT eae 
.Goul véeteh S20 iigneel,) 


2982 me, 
: Cnuge: OBGO ROGER, 


ee ee ee ee ee Sak a ell mi 


SHIR TURKS aA 


nga 
92540 G) aAy oF, rast 
Ceol: Dep. ceges Lay" vied 


LOC 

0129 
012A 
012B 
Or2ec 
012D 
012E 
012F 
0130 
0131 
0132 
Gtss 
0134 
ON35 
0136 
0137 
0138 
0139 
013A 
013B 
013C 
013D 
013E 
013F 
~—0140 


COSMAC CODE 
4F 
50 
59 
20 
2A 
4D 
2)! 
4F 
30 
OZ 
43 
45 
2A 
20 
2D 
44 
41 
54 
4 1 
28 
2A 
4C 
29 
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SOURCE CODE 


pete eO 
port Te. B 
(Rowe. oY 
nha ~«eSPACE 
pi cAee t ae 
fa. &M 
Hoes. «tS 
pha 6. $0 
BE Be) ay 
wee ©. SR 
pipes). 0. SG 
,#45 o os 
,#2A 0 oe 
ni 20-s —+2SPAGE 
SEZ > Bee. = 
pHoA- 5 8.D 
#41 A 
pe aed 
PPE SOA 
THOS & 1c ( 
(#OAS 2 o® 
pHaC 9G 
29.5 99.) 


| 
] 
| 


LOC COSMAC CODE 


8091 
8091 
8091 
8091 
8091 
8091 
8091 
8091 
8091 
8091 
8091 
8091 
8091 
8091 
8091 
8094 
3007 
8098 
8099 
809A 
809B 
SOC 
809E 
80A1 
80A3 
B0A4 
80A5 
80A6 
80A8 
80A9 
B0AA 
80AC 
80AD 
80AE 
80AF 
80B0 
80B1 
80B2 
80B3 
80B5 
80B6 
80B7 
80B8 
80B9 
80BA 
80BC 
80BD 
80BE 
80BF 
80C0 
80C2 


F86FAC 
F840B9 
93 

F6 

DC 

20 

99 
3A06 
F810A7 
F808A9 
46 

B7 
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SOURCE CODE 
- WRITE TO TAPE ROUTINE 
--FROM RCA VIP-711 OPERATING SYSTEM. 
REGISTER DESIGNATIONS: 
a 34°. Program counter 
wlen6 Pointer to next data byte 
ay 2 Parity register 
--R7.1 Holds byte being written 
eee Sync tone timer 
oeo.0 Battcounter 
eRe P counter to Bit-write subroutine 
--RE Page counter (last byte makes it 0) 


WRETB: 76E->RC.0 ..initializations, ~RC=4816F 

#40->R9.1 
SYNC:GHi 2 --otart’Synestone: i send out 

SHR ssarhighebat. 
SEP C 
DEC 9 ..-Decrement sync counter and 
GHI 9 .-continue sending high bits 
BNZ SYNCoT..fort5¥seconds. 

NEWBY:#10->R7.0 ..Reset parity register (Odd). 
#08->R9.0 
LDA R6 ..Get next data byte to send 
PHI 7 
GHI 3 feoena Out a Start bit 
SHL 
SEP C 
GLO 6 .-Is this the end of the page? 
BNZ NXTBT ..I£ yes, then decrement the 
DEC E wa page sCOUNTEr, 

NXTBT:GHI 7 ..-Get remaining data byte 
SHR ..and write next bit to tape. 
PHI 7 
SEP € 
DEC 9 MiSe CNL Seenemlast: Die: 
GLO 9 
BNZ NETEP?..liono;hgetinextebit. 
INC 7 .-If yes, finish up the 
GLO 7 .»parity calculation and write 
SHR «s Keatottapes 
SEP C eo 
GLO E .-Is this the last byte? 
BNZ NEWBY ..If no, then get new byte, 
SEP C ..else send a stop bit. 

VIDEO: , #69 ..Turn video display on. 
DEC 6 .-Point at last byte written. 
SEP 4 .-Call video interrupt routine. 

END1:BR END1 ..Terminate in this loop. 
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Taps 
on 
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Do. 


sacruon 3 
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LOC COSMAC CODE 


80C2 
80C2 
80C2 
80C2 
80C2 
80C2 
80C2 
80C2 
80C2 
80C2 
80C2 
80C2 
80C2 
80C5 
80C8 
80C9 
80CB 
80CC 
80CD 
80CF 
80D0 
80D2 
80D6 
80D7 
80D8 
80D9 
80DA 
80DB 
80DC 
80DE 
80DF 
80E0 
80E2 
80E3 
B80E4 
80E5 
80E6 
80E7 
80E9 
80EA 
80EB 
80ED 
80EF 


F883AC READ: 
F80AB9 AGAIN: 


3BCF 


TONE: 


WAITST: 


F809A9A7 


a7, 
76 
B7 
29 
DC 
89 
3AD6 
87 
F6 
33E3 
7B 
97 
56 
16 
86 
3ACF 
2E 
8E 
3ACF 
30BD 


NEXT: 


PARTY: 
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SOURCE LINE 


TAPE READ ROUTINE 
FROM RCA VIP-711 OPERATING SYSTEM. 


R6 Pointer to data buffer 


-R7.0 Partiy register 


R7.1 Temporarily stores incoming byte 


~R9 Sync tone timer 
-R9.0 Bit counter 


RC Read-a-bit P Counter 
RE.O Page counter 


#83->RC.0 ..Initialize, RC=#8183 
#OA->R9.1 


SEP C ».-Read ina bit 

BDF AGAIN ..Keep trying until a i-bit. 
DEC 9 .-Are there 2560 high bits 
GHI 9 se Nal row? 

BNZ TONE 

SEP C ..Yes, SO it must be a sync; 


BNF WAITST..Just wait for a start bit. 
#09->R9.0,R7.0 


GHi-7 ..Get data byte so far and 
SHRC meshiftwon theinext Dit 
PHites/ 

DEC 9 

SEP C .-Read in next bit. 

GLO 9 --Is this the last bit? 

BNZ NEXT ..No, then get the next one. 
GLO 7 t.Wes, ‘so check#ithe? parity. 
SHR 

BDF PARTY ..Turn Q on if there is 

,#7B »ea parity error. 

GHI 7 .-Store finished byte in 

STR 6 ..-the next buffer location 
INC 6 

GLO 6 .-Is this end of a page? 

BNZ WAITST ..No, then read next byte. 
DEC E »-ves, decrement page counter. 
GLO E »-Is this the very last byte? 


BNZ WAITST ..No, then get next else 
BR #VIDEO ..goto VIDEO and terminate. 


—. 
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LOC COSMAC CODE SOURCE LINE 


816E ie a a a es ee 
816E -- SUBROUTINE TO WRITE A BIT TO TAPE. 
2s .»-FROM RCA VIP-711 OPERATING SYSTEM. 
E SS ad et al a aaa atl 
816E - REGISTER DESIGNATIONS: 
816E --R3  P counter of calling program 
816E »-R7 Parity register 
816E --RC P counter of this subroutine 
816E sone timing counter 
816E Orn ere 
816E D3 RET1:SEP 3 .-Return to Tape-write. 
816F F80A LDIgeOA \.cltiitlalszertamer to.l0.. 
Si, 3B76 BNF ZERO ..Is this a zero bit? 
B73 F820 LDI #20 ..I1f£ yes, then let timer be 32 
Sabe 17 ENG 7 ..and increment parity reg. 
8176 7B ZERO: , #78 seeLUCNecOMOn eS CartyOote Dit. 
Sei © BE PHI F 
B178 FEO. MINUS:SMI #01 ..Count down timer with 
817A 3A78 BNZ MINUS..Q turned on. When timer 
817C 396E BNO RET1 ..1S zero, turn Q off and 
817E 7A ,#7A ..repeat timer loop with Q 
817F 9OF GHEE . OFF. When finished 
8180 3078 BR MINUS ..return to calling routine. 
8182 a ee ma a eg Sa ro 
8182 »-e SUBROUTINE TO READ IN A BIT FROM TAPE 
8182 »eFROM RCA VIP-711 OPERATING SYSTEM 
8182 PRE cath kant el oy eae a a oe ans aoe ian La a oar reba eT 
8182 ..REGISTER DESIGNATIONS: 
8182 Stet P counter of calling program 
8182 eR Parity register 
8182 RC jee acounter of jthis subroutine 
8182 ee eee eee 
Sio2 D3 REV Zt obeus .-Return to Tape-Read. 
8183 F810 EDL £10. Bawdnicialsze, Timer to -10;. 


8185 3D85 WAITSBN2 WAIT ..Wait for start of bit. 
8187 3D8F TIMER:BN2 OUT ..I£ the loop is exited from 


8189 .ethis point then it must 
8189 .-be a low bit. 

8189 FFO1 SMI #01 ..Continue timing. 

818B 3A87 BNZ TIMER..I£ we fall through the 
818D 17 ING 7 ~.timing loop it) must -be 
818E .-a high bit. Increment the 
818E »eparity register. 

818E 9C GHI iC ..-Get current byte. 

818F FE OUT: , #FE ..(SHL)Shift out this bit. 
8190 3590 WAIT2:B2 WAIT2 ..Wait for bit transition 
8192 3082 BR RET2 ..if necessary then return. 
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